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Abstract 

Canonical formulation of quantum field theory on the Light Front (LF) is re- 
viewed. The problem of constructing the LF Hamiltonian which gives the theory 
equivalent to original Lorentz and gauge invariant one is considered. We describe 
possible ways of solving this problem: (a) the limiting transition from the equal-time 
Hamiltonian in a fast moving Lorentz frame to LF Hamiltonian, (b) the direct com- 
parison of LF perturbation theory in coupling constant and usual Lorentz-covariant 
Feynman perturbation theory. The results of the application of method (b) to QED- 
1+1 and QCD-3+1 are given. Gauge invariant regularization of LF Hamiltonian 
via introducing a lattice in transverse coordinates and imposing periodic boundary 
conditions in LF coordinate x~ for gauge fields on the interval \x~\ < L is also 
considered. 
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1 Introduction 



A possible approach to solving the strong interaction field theory is the canonical quantization 
on the Light-Front (LF) with the application of corresponding Schroedinger equation. To realize 
such program one introduces LF coordinates [1] x ± = (x° ± x 3 )/^/2, x 1 , x 2 , where 
are Lorentz coordinates. The x + plays the role of time, and canonical quantization is carried out 
on a hypersurface x + = const. The advantage of this scheme is connected with the positivity of 
the momentum P_ (translation operator along x~ axis), which becomes quadratic in fields on 
the LF. As a consequence the lowest eigenstate of the operator P_ is both physical vacuum and 
the "mathematical" vacuum of perturbation theory [2]. Using Fock space [3] over this vacuum 
one can solve stationary Schroedinger equation with Hamiltonian P + (translation operator 
along x + axis) to find the spectrum of bound states. The problem of describing the physical 
vacuum, very complicated in usual formulation with Lorentz coordinates, does not appear here. 
Such approach is called LF Hamiltonian approach. It attracts attention for a long time as a 
possible mean for solving Quantum Field Theory problems. 

While giving essential advantages, the application of LF coordinates in Quantum Field 
Theory leads to some difficulties. The hyperplane x + = const is a characteristic surface for 
relativistic differential field equations. It is not evident without additional investigation that 
quantization on such hypersurface generates a theory equivalent to one quantized in the usual 
way in Lorentz coordinates [4-12]. This is in particular essential because of the special diver- 
gences at p_ = appearing in LF quantization scheme. Beside of usual ultraviolet regulariza- 
tion one has to apply special regularization of such divergences. We will consider the following 
simplest prescription of such regularization: 

(a) cutoff of momenta p_ 

\P-\ >£, £ > 0; (1) 

(b) cutoff of the x~ 

-L<x~<L. (2) 

with periodic boundary conditions in x~ for all fields. 

The regularization (b) discretizes the spectrum of the operator P_ (p_ = irn/L, where n is 
an integer). This formulation is called sometime "Discretized LF Quantization" [13]. Fourier 
components of fields, corresponding to p_ = (and usually called "zero modes") turn out to 
be dependent variables and must be expressed in terms of nonzero modes via solving constraint 
equations (constraints). These constraints are usually very complicated, and solving of them is 
a difficult problem. 

The prescriptions of regularization of divergences at p_ = described above are convenient 
for Hamiltonian approach, but both of them break Lorentz invariance and the prescription (a) 
breaks also the gauge invariance. Therefore the equivalence of LF and original Lorentz (and 
gauge) invariant formulation can be broken even in the limit of removed cutoff. To avoid this 
nonequivalence some modification of usual renormalization procedure may be necessary, see 
for example [11, 14] and [15]. The problem of constructing the LF Hamiltonian which gives a 
theory equivalent to original Lorentz and gauge invariant one turned out to be rather difficult. 
We will describe possible approaches to this problem. 
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In sect. 2 we give basic relations of quantum field theory in LF coordinates. In sect. 3 
we consider the limiting transition from fast moving Lorentz frame to the LF. This transi- 
tion relates Hamiltonian formulations in Lorentz and LF coordinates and firstly was presented 
in [7, 8]. In sect. 4 we investigate the relation between LF perturbation theory in coupling 
constant and usual Lorentz-covariant Feynman perturbation theory. With the help of this in- 
vestigation we show how to construct LF Hamiltonian giving a theory perturbatively equivalent 
to original one for Yukawa model, for QCD in four-dimensional space-time, and for QED in 
two-dimensional space-time (originally this was considered in [11, 14, 16-18]. In sect. 5 we 
carry out gauge invariant ultraviolet regularization of LF Hamiltonian introducing a lattice in 
transverse coordinates x x ,x 2 and using, instead of transverse components of gauge fields, com- 
plex matrices in color space as independent gauge variables on the links of the lattice [19, 20]. 
We find that LF canonical formalism for gauge theories with this regularization avoid usual 
most complicated constraints connecting zero and nonzero modes. However, the quantization 
leads to Lorentz-noninvariant results. 



2 Formal canonical quantization of Field Theory 

on the Light Front and the problem of bound states 

In order to find the bound state spectrum in some field theory quantized on the LF the following 
system of equations is usually solved: 

p+l*) = p;i*>, (3) 

P_\V) =P!_\V), (4) 
Pl|*) = 0, (5) 
where P± = {Pi, P 2 }. The mass of bound state is equal to 

m = \jlP'+P'- (6) 
It was taken into account that nonzero components of metric tensor in LF coordinates are 

9-+ = 9+- = 1, 9u = 922 = -1. (7) 

The operators P_, P± are quadratic in fields, and the solution of equations (4), (5) is trivial. The 
problem is in solving the Schroedinger equation (3). Physical vacuum \ fl) is lowest eigenstate 
of the operator P_, and must satisfy the equations: 

(8) 
(9) 
(10) 

To fulfil equations (8), (9) one should subtract , if it is necessary, corresponding renormalizing 
constants from the operators P + , P_. The plays simultaneously the role of mathematical 
vacuum of LF Fock space. A solution of Schrodinger equation (3) belongs to this space. 

Expressions for the operators P + , P_ can be obtained by canonical quantization on the LF. 
Let us describe the procedure of such quantization via some examples, without analyzing so far 
the question about the equivalence of appearing theory and original Lorentz-covariant one. It 
is assumed that in addition to explicit regularization of divergences at p_ = also ultraviolet 
regularization is implied. 



p + |n> 


= o, 


p_|n> 


= o, 


p ± \n) 


= 0. 
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2.1 Scalar selfinteracting field in (l+l)-dimensional space-time. 

Peculiarities of LF quantization are well seen even in this simple example. We have only LF 
coordinates x + , x~. The Lagrangian is equal to 

L — Jdx~ (d^ipd^ip - ^mV - V 4 ) , (11) 



or 



L = fdx~ (^d+tpd.tp - \m\ 2 - Ay? 4 ) , (12) 

The "time" derivative d + <p enters into this Lagrangian only linearly. For the transition to 
canonical theory it is sufficient to rewrite the expression / dx~ ^d + <pd-ip in standard canonical 
form. To achieve this let us take the Fourier decomposition 

oo 

<p(x~) = (2tt)~5 J dk \2k\~^ (a(k) exp(-ikx') + a + (k) exp(ifcaT)) , (13) 
o 

where k = (p(x~) = tp(x + ,x~), a(k) = a(x + , k). The Lagrangian (12) takes the form 



OO / 

V 



L = I dk \ ° + WKk)-a(k)a+(k)\ _ ^ (M) 



where a = da/dx + and 

H = jdx' QmV 2 + V) • (15) 
Here we have used the equality 

oo oo 

J dk J dk'5(k + k')k' (a(k)a(k') - a + (k)d + (k')) = 0. (16) 







It is implied that the function (p(x) in (15) is expressed in terms of a + (k) and a(k) with the help 
of formulae (12). "Time" derivatives a(k), a + (k) enter into Lagrangian L in a form standard 
for canonical theory. Therefore one can interpret after quantization the a + (k) and a(k) as 
creation and annihilation operators satisfying the following commutation relations at fixed x + 
and k > 0, k! > 0: 

[a(k), a + (k')} = 5{k - k'), [a(k), a{k')} = 0, (17) 

It is also seen that the H is LF Hamiltonian, i.e. H = P + . 
We have also the formulae 

P„ = J dx'T.^ (18) 
where the energy-momentum tensor T Vfl is equal to 

T vll = dvtpdptp - g^C. (19) 
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Via this relation one can reproduce the expression (15) for P + = H , and obtain the equality 

oo 

P- = J dx- (<9_y?) 2 = ^Jdkk (a + (k)a(k) + a(k)a + (k)) . (20) 

o 

The lowest eigenstate of the operator P_ is the physical vacuum \Q) for which 

a(fc)|n> = (21) 

at any k. It is seen that vacuum expectation values (Q|P_|fi), (Q|P+|Q) are infinite. The renor- 
malization can be got by taking normal ordered forms : P+ :, : P : with respect to the operators 
a + , a (the symbol :: means as usual that operators a + stand everywhere before of operators 
a). Normal ordering of the Ay? 4 term in the Hamiltonian P + leads also to renormalization of 
the mass. In the following, writing P + , P_, we mean the expressions : P + :, : P_ :, satisfying 
conditions (8), (9). Normal ordering of the operators P + , P_ allows to avoid all ultraviolet 
divergences in this simple model. 



2.2 Theory of interacting scalar and fermion fields 

in (3+l)-dimensional space-time (Yukawa model) 

The Lagrangian of the model is 

L = J d 2 x ± dx~ (i) (ry^ — M) ip + ^d„(pd^<p - V - ffW<P ~ A> 3 - V 4 ) , (22) 

where M is the fermion mass, m is the boson mass, ip = t/> + 7°, tp = <p + ; g, A, A' are coupling 
constants. Here and so on /z, z/, . . . = +,—,1,2; i,k,... = 1,2; x L = (x 1 ^ 2 ). For Dirac's 
7-matrices we use: 

o ( -il \ 3 / il \ ± ( -ia L \ , _. 

^={u o J' -?={ii o J' -r L= { ia- J' ^ 

where / is a unit 2x2 matrix, a 1 - = {a 1 , a 2 }, a 1 are Pauli matrices. 
We introduce 2-component spinors Xi writing 

*=({)■ (24) 

The Lagrangian L can be written in the form 

L = j d 2 x ± dx~ (iV2 X + d +X + i>/2f + (a' l di -M)x + h.c.) + 

+d + <pd-<p - - \m 2 p 2 - igtp ~ X + t) - A> 3 - Ay? 4 ) , (25) 

where h.c. means Hermitian conjugation. The variation of this Lagrangian with respect to x + 
leads to the equation 

V2d^ = - (a% -M)x + gpx- (26) 
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This equation does not contain derivatives in x + and therefore is a constraint. One should solve 
it with respect to £ and substitute the result into the Lagrangian. In doing this we must invert 
the operator d- which becomes an operator of multiplication ik- after Fourier transformation. 
Inverse operator [ik)~ x has singularity at = 0. To avoid this singularity we introduce the 
regularization (1). For any function f(x~) = f(x + ,x~,x x ) we define Fourier transform 



/OO 

where /(&-) = f(x + , /c^ar 1 ), and put 

1 



dk- e 



ik-x 



[/ooi = 



dk_e tk - x ~f{k_), |A;„|>£>0. 



(27) 



(28) 



We insert into the Lagrangian (25) the variables [%], [£], [£ + ], [<£>] instead of x, X + > 
£ + , (/? and obtain the constraint equation 

(29) 



(30) 



y/2d-[^] = -(a i d i -M) [x]+g[ 
instead of (26). Its solution is 

[Z] = -j=dz 1 (-(* i a i -M) [x] + *][*]]), 

where the operator dz 1 is completely defined by the condition 

d- l [f] = Id-Jif] . 



(31) 



After Fourier transformation the operator dz 1 is replaced by {ik-)' 1 . 

Substituting the expression (30) into the Lagrangian (where all fields x, x + , • • • are replaced 
with [x], [x + ], • • •) we come to the result: 

L = J d 2 x L dx- (iy/2 [ X + ] d + [ X \ + d-[ip\d+[<p\ + 

+-^( (°% ~ M) ix] - gM\x}}) + (-id-T l (((T k d k - M) ix] - g[[<p]\x]\) ~ 

-\dMm - \m 2 M 2 - A'M 3 - AM 4 ), (32) 

As in sect. 2a time derivatives 9+[x], d + [ip] enter into the Lagrangian (32) linearly. Therefore 
to go to canonical formalism it is sufficient to find a standard canonical form for the expression 

iV2\x + ]d + [ X } + d-[p + ]d + [p] (33) 
(before a quantization the quantities x + 5 X are elements of Grassman algebra). We write 

oo 

[<p(x-)\ = (27T)- 1 / 2 J dk_ (2k„y 1/2 (a(k„) exp(-iJfe_aT) + a+(A;_) exp(ik-x~)), (34) 

£ 

OO 

[Xr{x~)\ = (2tt)- 1/2 2- 1/4 J dk-(b r (k-) exp(-ik-x~) + c+(£;_) exp(ik-x~)), (35) 
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[x+(aT)] = (2tt)- 1 / 2 2^ 1 / 4 J dk-(c r (k-) exp(-ik-x-) + b+(k-) exp(ik-x-)), (36) 

£ 

where 

[</?(a; _ )] = [ip{x + , x~ , or 1 )], a(/c_) = a(x + , x -1 ) (37) 
et cetera, r = 1,2. The Lagrangian (32) takes the form 

oo 

d 2 x ± J dk^ (a(k^)d + (k^) - a + (k_)d(k_) - 

£ 

-b r (k-)b+(k-) - b?(h-)b r (k-) - c r (A;_)c+(A;_) - c+(k_)c r (k^)) - H, (38) 

where 

H = ( - ^( (a% - M) [x] - 5[M[x]]) + (-^_)- 1 ( (a k d k - m) [ X ] - <?[[*]]) + 

+\dMm + \m 2 [pf + A^f + A^ 4 ). (39) 

It is assumed that the quantities [x], [% + ], [<£>] m the formulae (39) are expressed in terms of 6, 
6+, c, c+, a, a+ with the help of (34), (35), (36). 

It follows from (38) that a + , a, b + , b, c + , c play a role of creation and annihilation operators. 
After quantization they satisfy the commutation relations (at x + =const): 

[a(fc_, ar 1 ), a + (A;'„, x ,J -)]_ = <S(*- - k'J)5 2 {x L - x f± ), (40) 
[6(fc_, x x ), 6 + (A; , _, x ;± )] + = 5(A;_ - k'_)5 2 {x ± - x f± ), (41) 
[c(fc_, or 1 ), c + (k'_, x /J -)] + = <S(Jfe_ - ^)5 2 (x x - x' x ), (42) 

where [x, y]± = xy ± yx. The remaining (anti) commutators are equal to zero. It is seen that 
the quantity H is LF Hamiltonian (H = P+). The operator of the momentum P_ is equal to 

oo 

P- = j d 2 xdx~ T__ = -J d 2 x ± J dk_ k4a + (k_)a(k„) + a(k_)a + (k_) + 

£ 

+b+(k_)b r (k_) - b r (k_)b+(k_) + c+(k_)c r (k_) - c r (k_)4(k_)). (43) 

The quantities P + = H and P_ should be normally ordered with respect to creation and 
annihilation operators. The lowest eigenstate of the momentum P_ is the physical vacuum. It 
is defined by conditions 

a(k-,x 1 -)\Sl)=Q, b(k-,x ± )\Q) — 0, c(Jfe_, x ± )|Q) = 0, Vi 1 ,L>e>0. (44) 

The equalities (8), (9) become true after normal ordering of the operators P + and P_. For the 
P_ it is seen from the formulae (43), and for the P + it follows from the fact that every term 
of P + contains a 5-function of difference between the sum of momenta k- of creation operators 
and the sum of momenta fc_ of annihilation operators. Due to the positivity of all momenta k-, 
in our regularization scheme every term of the P + contains at least one annihilation operator. 
Therefore for normally ordered P + we have P+|f2) = 0. 

The model under consideration requires ultraviolet regularization. It can be done in different 
ways. We consider this question together with the renormalization problem in sect. 3 and 4. 
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2.3 The U(N)-theory of pure gauge fields 

We consider the U (N) rather than the SU(N) theory because it is technically more simple. The 
transition to the SU (N) can be done easily. Gauge field is described by Hermitian matrices 

A,(x) = A;(x) = {A*(x)} , (45) 

where /! — +,—, 1, 2; i,j = 1,2, ... , N. Let us assume, that for the indexes i,j and analogous 
the usual rule of summation on repeated indexes is not used, and where it is necessary the sign 
of a sum is indicated. Field strengths tensor is 

F/iu = d^A v - d u A fl - ig[A», A v ], (46) 

and gauge transformation has the form 

A, - 4 = U+A^U + -U + d,U, U + U = I. (47) 

To escape a breakdown of gauge invariance we apply the regularization of the type (2) with 
periodic boundary conditions 

A li (x + ,-L,x ± )=A li (x + ,L,x ± ), (48) 

on the interval —L < x~ < L (this exact periodicity can be always achieved starting from 
gauge invariant one, i. e. the periodicity up to a gauge transformation). All Fourier modes of 
A^x) in x~ must be kept, including zero modes (at A;_ = 0). 
The Lagrangian has the form 

L = ~2~J d2x± I dx ' Tl ( F ^ F ^ ' ( 49 ) 



-L 



or 



L = J d 2 x^ J dx-Ti [f 2 + _ + 2F_ k F +k - ^F kk ,F kk ^ , (50) 

— L 

where k, k' = 1, 2. Time derivatives d + A k are present only in the term 

Tr (2(d„A k - d k A_ - tg[A_, A + ])d + A k ) . (51) 

This expression can be rewritten in standard canonical form only after fixing the gauge in a 
special form of the type A_ = because then the term above becomes similar to that of scalar 
field theory Tr (2(d^A k )d + A k ). However not every field, periodic in x~, can be transformed to 
the A_ = gauge. Indeed, the loop integral 

f ( f \\ 

T(x + ,x ± ) = Tr I Pexp % \ dx~ A_(x + , x~ , x x ) \\ , (52) 

V ~L J) 
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where the symbol 'P' means the ordering of operators along the x~, is a gauge invariant quantity. 
If this integral is not equal to N for some field then the gauge A_ = is not possible for this 
field. Therefore we choose more weak gauge condition 



and put 



A i3 _ = 0, i± j, d_A* = 0, 



A il _{x) = v i (x + ,x ± ). 



(53) 



(54) 



The gauge (53) breaks not only the local gauge invariance but, unlike the A_ = gauge, also 
global gauge invariance (it remains only the abelian subgroup of gauge transformations not 
depending on x~). This has some technical inconvenience but now any periodic field can be 
described in the gauge (53). Furthermore, if we restrict the class of possible periodic fields by 
the condition A_ = 0, disregarding the described consideration, we come to canonical theory 
with even more complicate constraints if zero modes are taken into account [21, 22, 24]. 

From the point of view of LF canonical formalism the variables A % 1 are "coordinates". 
Therefore one can restrict their values by the condition (53) directly in the Lagrangian without 
loosing any equations of motion. Let us introduce the denotations 



D-A% = (cL - ig{v l - v j ))A%, 



(55) 



where obviously 



D-A i > = {d--ig{v i -vi))At 



Also for any function f(x ) = f(x + ,x ,ar L ) periodic in x we denote 



(56) 



(57) 



i r 

/(o) = ^ J dx f(x ), 



(58) 



[/(*")] = f(x-) - / (0) . 



Obviously, 



(59) 



dx~ [f(x-)\ = 0. 



(60) 



After gauge fixing (53) the Lagrangian (50) takes the form 



L = J d 2 x J dx-hj2(d- 

t \ i 



d k d- 



Al 



Alt 



+ 2 £ (p_Af) d+Aj! + 2 (d- [Al\f + 



i,3, i^j 



~ i9 E (^'D-AC - (D_Af) AC) 



+ 
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+2 £ ^ LlLAf - ig E (4^-4" " (^-40 4 1 )) ~\Z F^F^\-\ 

L 

+^ / dx- E 0^-4* - 0-40 4 



(61) 



Here we have ignored some unessential surface terms. 

Variation of the Lagrangian in [A™] , A\ at i ^ j leads to constraints, the solution of which 
can be written in the form 



[A%] = d 



-2 



d k d-[At] -ig E (4'd-aC - (D_Af)A£) 

j',j'¥=i 



A l { 



i+3 



DZ 2 (d k D_A« - igY,{4D-4 j ~ (D-Af)A^) 



(62) 



(63) 



The operator d_} is completely defined, as before, by the condition (31) being well defined 
on functions [/(x)]. The operator Dz 1 after Fourier transformation in x~ is reduced to the 
multiplication by (i(k- — g{v l — v J ))) 1 . Therefore it has, in general, no singularities for any 
k- = n(n/L) with integer n. 

Substituting the expressions (62), (63) into the Lagrangian (61), we exclude from it the 
quantities [A™] and A % 1 at i ^ j. The variation of the Lagrangian in ^4+( ) leads to the 
constraints 

/ k \ 



Q ii (x + ,x ± ) = -2 



2Ld k d k v* + J dx- E {4'D-4-(D-4)4 i ) 



0. (64) 



They are first class constraints which can be posed on physical state vectors after quantization. 
Therefore we can keep the term with this constraint in the Lagrangian. 

Now we must put in the standard canonical form the terms of the Lagrangian 

L 



j d\ J dx- 2E(d-[4i)d + [4i + 2 E (D-A^d+A 



i,3, i+3 



It can be reached by going to Fourier transform 
1 



[A*(x-)] = —= E kZ 1/2 {al(k_)eM-ik-X~)+al + (k_)eMik-X-)}, 



2V2L k _ =n/L 
where we sum over A;_ = nn/L, n = 1,2,..., and 



(65) 



(66) 



¥i 2 



4^ E (k- - g{v l - v 3 )) 1/2 a k ] + (k^)exp(ik^x ) + 

v 2L {k->g(vi-vi) 



1/2 

+ E ( k - ~ 9( yj - yi )) a£(k-)exp(-ik-x~) 

k->g(vi — v l ) 



(67) 
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where we sum over all k_ = nn/L satisfying corresponding inequalities. The expression (65) 
takes the form 



/( ^ oo 
d 2 ^ E E (al(k„)d + (al) + (k„)-al + (k_)d + al(k_)) + 

{ i fc_=7r/L 

+ E E {4{k-)d + 4\k_)-4\k_)d + 4{k_))\. (es) 

Further, in the Lagrangian (61) there is a part 

L V = J rfV (2LE0VT -4LE(^(0))9 + ^} • 



The "momentum" conjugated to t>* is 



Hence, 



d + v l = —V* + d k A% 0) . 



The corresponding part of the Hamiltonian equals to 

H v = j E {V'd+v') -L V = J £x^2L E + ^(0)) 
and the corresponding part of canonical Lagrangian is 

Lv = y^E^^) -Hv 



(69) 



(70) 



(71) 



(72) 



(73) 



Excluding from the Lagrangian (61) the quantities and A % 1_ (at i 7^ j) via the equations 
(66), (67), replacing the terms (65) by the expression (68) and the part (69) by the expression 
(73), we obtain the result 



d 2 ^ E E (4(M<9 + 4 + (A:- 

[ i k-=7T/L 



) - af(M<9 + 4(M + 



+ E E (4(k-)d + 4 + (*_) - aSf + (*-)9 + aS?'(*-)) + 

+ ^ + £ ^ (0) q4 - If, (74) 

i i J 

where Q M are defined by (64), and the Hamiltonian H = P + is equal to 

H = fdFx fdx-lz(d-[A«]) 2 + E (D_A^D_A^- 1 -J2F k ^F kl A + 

+2Lf^x L '£(^r + d k At 



AL 



l fc(0) 



(75) 
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It is implied that instead of the quantities [A™] and A % 1 (at % + j) one uses the expressions 

(62), (63) and the A k j are expressed in terms of a\, a\ + , a% , a k j+ , (at i + j) and of A k ^ with 
the help of equations (66), (67) and 



Al 



Al 



+ A 



fc(0)- 



(76) 



It is seen from the formulae (74) that a l k + , a k , a% , play the role of creation and 
annihilation operators. After quantization they satisfy the following commutation relations (at 

x + = const): 

4(fc_,^),af (fc^)! =^5 kl 5 k _ >k ,_5 2 (x ± -x' ± ), (77) 



Also we have 

V i (x ± ),v j (x' ± )]_ = -i5 ij 5 2 (x ± - x' x ) 

Remaining commutators are equal to zero. 

The operator of the momentum P_, defined by 



8 ii '6tf6 kl 6 k _, k >_6 2 (x ± -x ,± ), i + j, i'+f- 



(78) 



(79) 



P- = j d 2 x j dx'T—, 



(80) 



acts on physical states |\&), satisfying the condition 

Q u ( X ± m = o, 

equivalent to the canonical operator 

=Jd 2 x\Y J E k.a k + (k_)al(k.) + J2 E 

\ i k^=w/L hj,i^j k->g(v i 

where the normal ordering was made. 

Physical vacuum |Q) satisfies the relations 

al(k.,x ± )\n) = 0, 



(81) 



fc_ar(*>!?'(*-) , (82) 



(83) 



al(k.,x ± )\n) = 0, i + j, 



(84) 



and the condition (81). 

This scheme is connected with the following essential difficulty The zero modes A^Jx 1 ) 
are present in the Lagrangian (74) and in the Hamiltonian (75) but the derivatives d + A k ^ are 
absent there. Therefore new constraints arise 



5H 

^ ( o)(^) 



= 0. 



(85) 
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These constraints are of the 2nd class and they must be solved with respect to ^(o) an d then 
the A % y Q \ have to be excluded from the Hamiltonian. The constraints (85) are very complicated 
and explicit resolution of them is practically impossible. The application of Dirac brackets does 
not simplify this. 

Due to this difficulty a practical calculation usually ignores all zero modes from the be- 
ginning. It makes the approximation worse. It is interesting that in the framework of lattice 
regularization it is possible to overcome the difficulties caused by the constraints (85) [19, 20]. 
This question will be considered in sect. 5. 

3 Limiting transition from the theory in Lorentz 
coordinates to the theory on the Light Front 

To clarify the connection between the theory in Lorentz coordinates in Hamiltonian form and 
analogous theory on the LF we perform the limiting transition from one to the other. Here this 
transition is considered in the fixed frame of Lorentz coordinates by introducing states that 
move at a speed close to the speed of light in the direction of the x 3 axis. Constructing the 
matrix elements of the Hamiltonian between such states and studying the limiting transition to 
the speed of light (an infinite momentum), we can derive information about the Hamiltonian 
in the light-like coordinates. This information also takes into account the contribution from 
intermediate states with finite momenta. Here, we illustrate the results of such an investigation 
using (l + l)-dimensional theory of scalar field with the A<^ 4 -interaction. Instead of x 3 we denote 
analogous space coordinate by x 1 . The generalization of the method to (3 + l)-dimensional 
Yukawa model is discussed briefly at the end of this section. The limiting transition studied 
here is accomplished approximately by subjecting the momenta p\ to an auxiliary cutoff that 
separates fast modes of the fields (with high pi values) from slow modes (with finite p\ values). 
This cutoff is parameterized in terms of the quantities A, Ai, and 5 and the limiting-transition 
parameter r] (rj > 0, r] — > 0): we have ?? -1 Ai > \pi\ > r]~ 1 5 for the fast modes and p\ < A 
for the slow modes (A ^> 5). For i] 0, the inequality r]~ 1 5 > A holds, so that the above 
momentum intervals are separated. The field modes with the momenta i]~ 1 5 > \p±\ > A 
are discarded. This procedure is justified by the fact that the resulting Hamiltonian in the 
limit ?7 — > reproduces the canonical LF Hamiltonian (without zero modes) when only the 
fast modes are taken into account and is consistent with conventional Feynman perturbation 
theory for 5^0. Therefore, even an approximate inclusion of the other (slow) modes may 
provide a description of nonperturbative effects, such as vacuum condensates. The effective 
LF Hamiltonian obtained here for the model under consideration differs from the canonical 
Hamiltonian only by the presence of the vacuum expectation value of the scalar field and by 
an additional renormalization of the mass of this field. The renormalized mass involves the 
vacuum expectation value of the squared slow part of the field. Masses of bound states can be 
found by solving Schrodinger equation 

2 

Tfl 

P+\*) = 2^l*>, ( 86 ) 

with obtained Hamiltonian P + . 

We start from the standard expression for the Hamiltonian of scalar field (p(x) in (1 + 1)- 
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dimensional space-time in Lorentz coordinates x^ = (x ,^ 1 ), at x° = 0: 



H 



2 (^) 2 + ^ 2 + Av^ 4 



(87) 



where Il(a; 1 ) are the variables that are canonically conjugate to tp{x l ) = tp(x° = O^ 1 ), and 
the symbol : : of the normal ordering refers to the creation and annihilation operators a and 
a + that diagonalize the free part of the Hamiltonian in the Fock space over the corresponding 
vacuum |0). These operators are given by 



^(x 1 ) = — / dpi(m 2 + p\)~ 1 ^ a(p\) exp(— ipix 1 ) + h.c. 
J An J 1 



(88) 



n(x x ) = -=L f d Pl {m 2 + pl)~ 1/4 fa(pi) exp(-ip 1 a; 1 ) - h.c] , (89) 



47T 

where a(pi)|0) = 0. 

To investigate the limiting transition to the LF Hamiltonian (defined at x + = 0), it is more 
convenient to go over from the Hamiltonian (87) to the operator H + Pi = \^2P + , where the 
momentum Pi has the form 

Pi = J dpia + {pi)a{pi)pi. (90) 

Applying the above parametrization of high momenta in terms of rj, i] — > 0, to p\ we can then 
consider the transition to an infinitely high momentum of states as a limit of the corresponding 
Lorentz transformation with parameter r/. To be more specific, we have pi — > (—i]^/2)~ 1 q^, 
where g_ is a finite momentum in the light-like coordinates, and 

Um^v^)- 1 ^!^ + Pi) a! o = obi)) = (q'_\(P + ) x+=0 \q-). (91) 

It follows that the eigenvalues E + of the operator (P + ) x + =0 that correspond to the momentum 
g_ are obtained as the corresponding limit of the eigenvalues E{rj) of the operator (H + Pi) x o =0 
at momentum pi. 

E + = lim(77v / 2) _1 J B(?7). (92) 

In the following, we consider this limiting transition as a part of the eigenvalue problem for 
the operator H + Pi, using perturbation theory in the parameter rj. Separating the Fourier 
modes of the field into fast and slow ones, as is indicated above, and neglecting the region of 
intermediate momenta (?7~ 1 <5 < \pi\ < i]^ 1 Ai is the region of the fast modes, and \pi\ < A is 
the region of the slow modes), we can substantially simplify this perturbation theory. The 77 
dependence of the field operators and Hamiltonian can then be determined by making, in the 
region of fast momenta, the change of the variables as 

Vl = i 1 - 1 k, a(pi) = y/rja(k), 5 < \k\ < A x , [a(k), a + (k')\ — 8(k - k'). (93) 

The fast part [ip(x 1 )]f of the operator ^(x 1 ) is estimated as 

[<p(x 1 )] f = <p(y) + 0(r j 2 ), y = r j - 1 x 1 , (94) 
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<p( y ) = (47T)- 1 / 2 / dk\k\~ 1/2 [a(k) exp(-iky) + h.c.]. (95) 



We denote the slow part of the field f by <p (if = [ip\f + <p). Substituting formulas (93)-(95) 
into Hamiltonian (87), we obtain 

H + P 1 = r]- 1 (h + 77/n + rfh 2 + ...), (96) 



"1 

h = 2j dkd + (k)d(k)k, (97) 



h 1 = (H + P 1 )^ n=n = (H + P 1 ), (98) 



h 2 = J dk ( ^ j ~a + (k)~a(k)+ : A J dy [<p 4 (y) + 4£(0)£ 3 (y) + 6y5 2 (0)^%) 



(99) 



<5<|fc|<A 



Prior to performing integration with respect to y, we formally expanded the operators ^(x 1 ) = 
<f(rjy) in Taylor series in the variable rjy and estimated their orders in the parameter 77 at fixed 
y. Such estimates can be justified at least in Feynman perturbation theory. The operator 
(H + Pi) in (98) is defined in such a way that its minimum eigenvalue is zero. Let us consider 
perturbation theory in the parameter rj for the equation 

(H + Pi)\f) = E\f) (100) 

under the condition that the states |/) have, as in formula (91), a negative value of Pi propor- 
tional to i]^ 1 for 77 — > and also describe the states with a finite mass. The expansions of the 
quantity E and the vector |/) in power series in 77 can then be written as 

00 00 

E = V ^J2v n En, \f) = T,V n \fn). (101) 

n=2 n=0 

We arrive at the system of equations 

Wo> = 0, (102) 



h o \fi) + h 1 \f ) = 0, (103) 

W2) + hi\h) + (h 2 - E 2 )\f ) =0, .... (104) 

To describe solutions to these equations, we use the basis generated by the fast-field operators 
d + (k) over the vacuum |0) and the slow-field operators <p and n over the vacuum \v) that corre- 
sponds to the Hamiltonian (H + Pi). The vectors of this basis can be symbolically represented 

as 

d + ...d + \0)fi...fifl...fl\v). (105) 
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By virtue of (97), the manifold of solutions |/ ) to equation (102) is reduced to the set of vectors 
(105), which do not contain the operators a + (k) with k > 5. Let Vq be the projection operator 
onto this set. According to equation (103), we then have 

V o ho\fi) = (H + Pi)\f ) = 0. (106) 

Equation (106) requires that the vectors |/ ) be the lowest eigenstates of the operator (H + P^, 
that is, linear combinations of basis vectors (105) including neither the operators a + {k) with 
k > 5 nor the operators <p and II. We denote the projection operator on this set of vectors (105) 
by Vq. To determine the quantity Ei which we are interested in, it is sufficient to consider the 
"Pp-projection of equation (104). Taking into account (102), (105), and (106), we find that E 2 
appears as a solution to the eigenvalue problem 

V h\f )=E 2 \f ). (107) 

Thus, in accordance with (91) and (92), the operator V^Yl^Pq pl a Y s the role of the effective LF 
Hamiltonian P+. Substituting formula (99) for the operator /i 2 into the expression forP^, we 
take into account that, between the projection operators Vq, the contribution of the field modes 
with positive momenta (k > 5) vanishes and that the products of the operators of the slow part 
of the field can be replaced with their expectation values for the vacuum \v). In addition, we 
note that, under the Lorentz transformation corresponding to the limiting transition n — > in 
formula (91), the variable y goes over into the light-like coordinate y~ = —y/ y/2, the momenta k 
go over into the light-like momenta g_ = —\^2k, and the corresponding coordinate y + vanishes 
at x° = (for finite values of y~). Going over to the operators A(q_) = 2~ 1 / 4 a(— k) for k < —5 
(<?- > 8\/2) and to the corresponding field $(y + = 0, y~) = (f>-{y), where y5_ is the part of 
the (p containing only the modes with negative momenta (k < —5), we obtain the effective 
Hamiltonian in the form 

$ 2 + 4A(y3)„$ 3 + A$ 4 j :, (108) 

where (. . .)„ is the expectation value for the vacuum \v). Expression (108) coincides with the 
canonical effective Hamiltonian if, in the latter, we take into account the shift of the field by 
the constant (<p) v and the change in the mass squared by 12A [(: 2 :) v — (0)1]. 

Analogous results were obtained for Yukawa model in (3 + l)-dimensional space-time [8] in 
regularization of Pauli-Villars type, introducing a number of nonphysical fields with very large 
masses. The absence of essential difference between the Hamiltonian obtained via limiting 
transition and canonical LF Hamiltonian is connected with this choice of regularization. Other 
regularizations can lead to more complicated results. 

This method of limiting transition can not be directly expanded to gauge theories, because 
the approximations used for nongauge theories are not justified. 

4 Comparison of Light Front perturbation theory 
with the theory in Lorentz coordinates 

As is already known, canonical quantization in LF, i.e., on the x + = const hypersurface, can 
result in a theory not quite equivalent to the Lorentz- invariant theory (i.e., to the standard Feyn- 
man formalism). This is due, first of all, to strong singularities at zero values of the "light-like" 



*ff 



dy 



m 



+ 12A(: 2 :) v 
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momentum variables Q_ = -^(Qo — Qz)- To restore the equivalence with a Lorentz-covariant 
theory, one has to add unusual counter-terms to the formal canonical Hamiltonian for the LF, 
H = P + . These counter-terms can be found by comparing the perturbation theory based on the 
canonical LF formalism with Lorentz-covariant perturbation theory [11]. This is done in the 
present section. The LF Hamiltonian thus obtained can then be used in nonperturbative calcu- 
lations. It is possible, however, that perturbation theory does not provide all of the necessary 
additions to the canonical Hamiltonian, as some of these additions can be nonperturbative. In 
spite of this, it seems necessary to examine this problem within the framework of perturbation 
theory first. 

For practical purposes a stationary noncovariant LF perturbation theory, which is similar to 
the one applied in nonrelativistic quantum mechanics, is widely used. It was found [25-27] that 
the "light-front" Dyson formalism allows this theory to be transformed into an equivalent LF 
Feynman theory (under an appropriate regularization). Then, by re-summing the integrands 
of the Feynman integrals, one can recast their form so that they become the same as in the 
Lorentz-covariant theory. (This is not the case for diagrams without external lines, which we 
do not consider here.) Then, the difference between the LF and Lorentz-covariant approaches 
that persists is only due to the different regularizations and different methods of calculating 
the Feynman integrals (which is important because of the possible absence of their absolute 
convergence in pseudo- Euclidean space). In the present section, we concentrate on the analysis 
of this difference. 

A LF theory needs not only the standard UV regularization, but also a special regularization 
of the singularities Q- = 0. In our approach, this regularization (by method (1)) eliminates 
the creation operators a + (Q) and annihilation operators a(Q) with \Q l _\ < e from the Fourier 
expansion of the field operators in the field representation. As a result, the integration w.r.t. 
the corresponding momentum Q_ over the range (— oo, —e) U (e, oo) is associated with each line 
before removing the (^-functions. Different propagators are regularized independently, which 
allows the described re-arrangement of the perturbation theory series. On the other hand, this 
regularization is convenient for further nonperturbative numerical calculations with the LF 
Hamiltonian, to which the necessary counter-terms are added (the "effective" Hamiltonian). 
We require that this Hamiltonian generate a theory equivalent to the Lorentz-covariant theory 
when the regularization is removed. Note that Lorentz-invariant methods of regularization (e.g., 
Pauli-Villars regularization) are far less convenient for numerical calculations and we shall only 
briefly mention them. 

The specific properties of the LF Feynman formalism manifest themselves only in the inte- 
gration over the variables Q± = -^(Qo ± Q3), while integration over the transverse momenta 
Q ± = {Qi,Q2} is the same in the LF and the Lorentz coordinates (though it might be non- 
trivial because it requires regularization and renormalization). Therefore, we concentrate on a 
comparison of diagrams for fixed transverse momenta (which is equivalent to a two-dimensional 
problem). 

In this section we propose a method that allows one to find the difference (in the limit 
e — > 0) between any LF Feynman integral and the corresponding Lorentz-covariant integral 
without having to calculate them completely. Based on this method, a procedure is elaborated 
for constructing an effective LF Hamiltonian correct to any order of perturbation theory. This 
procedure can be applied to nongauge field theories as well as to Abelian and non-Abelian 
gauge theories in the gauge = with the gauge vector field propagator chosen according 
to the Mandelstam-Leibbrandt prescription [28, 29]. The question of whether the additions 
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to the Hamiltonian, that arise under that procedure, can be combined into a finite number of 
counter-terms must be considered separately in each particular case. 

We will consider the application of this formalism to Yukawa model, to QCD in four- 
dimensional space-time and to QED in two-dimensional space-time. 

4.1 Reduction of Light Front and Lorentz-covariant Feynman 
integrals to a form convenient for comparison 

Let us examine an arbitrary one particle irreducible Feynman diagram. We fix all external 
momenta and all transverse momenta of integration, and integrate only over Q + and Q-\ 

F=lim/ n - rf2Q ' m '' pt) (109) 

We assume that all vertices are polynomial and that the propagator has the form 

ZiQ) « , m ; Z[Q) ± r, ■ , ("<>) 



Q 2 - m? + ias ' (Q 2 - m? + iae) (2Q+Q_ + iae) 

where z(Q) is a polynomial. A propagator of the second type in (110) arises in gauge theories 
in the gauge = if the Mandelstam-Leibbrandt formalism [28, 29] with the vector field 
propagator 

l ( QjX + QJl on \ 

9^ ~ ^ ^ — ~-^ 2 Q+ , (in) 



Q 2 + ias V ^ 2Q+Q_ + iae 

is used. In equation (109) either Mf = mf + Q^ 2 ^ 0, where m; is the particle mass, or Mf = 0. 

The function / involves the numerators of all propagators and all vertices with the necessary 
^-functions, that include the external momenta p k (the same expression without the 5-functions 
is a polynomial, which we denote by /). We assume for the diagram F and for all of its 
subdiagrams that the conditions 

w„ < 0, u + < 0, (112) 

hold, where oj + is the index of divergence w.r.t. Q+ at Q l _ ^ Vi, and cu n is the index of 
divergence in Q + and Q_ (simultaneously); Q± = -^(Qo ± Qz)- The diagrams that do not 
meet these conditions should be examined separately for each particular theory (their number 
is usually finite). We seek the difference between the value of integral (109) obtained by the 
Lorentz-covariant calculation and its value calculated in LF coordinates (LF calculation). 
In the LF calculation, one introduces and then removes the LF cutoff \Q- \ > e > 0: 

^ = & Ss / n **- / n «*> nim f\% + fa) ■ < 113 > 

Ve 

where 14 = llj(( — oo, — e) U (e, oo)). Here (and in the diagram configurations to be defined 
below) we take the limit w.r.t. e, but, generally speaking, this limit may not exist. In this case, 
we assume that we do not take the limit, but take the sum of all nonpositive power terms of 
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the Laurent series in e at the zero point. If conditions (112) are satisfied, Statement 2 from 
Appendix I can be used. This results in the equality 

F lf = lim lim fUdg* [ TT^" ^/if'^L ( 114 ) 

From here on, the momenta of the lines are assumed to be expressed in terms of the loop 
momenta q k , is a sphere of a radius L in the g^-space, and L depends on the external 
momenta. Now, using Statement 2 from Appendix I, we obtain 

r r f(O i n s ) p-fT ll Q i + 2 -i 3 Y ll Q- 2 

F u = lim lim lim lim / TT dq k / TT <k- Z*\ n „ Vv iTT? r^— . (115) 

^o^o^o 7 ^o y 11 y+ y 11 y " f]i(2(5 l + (5!_ - + iae) v ; 

To reduce the covariant Feynman integral to a form similar to (114), we introduce a quantity 
F: 

F = hm lim lim / TW /W ' P i ' 2 ^r— ■ (116) 

^0^07-oy 11 y n .(2Q* + Q*_ - M? + iae) v ; 

Let us prove that this quantity coincides with the result of the Lorentz-covariant calculation 
F cov . To this end, we introduce the a-representation in the Minkowski space of the propagator 



z{Q l ) 



_ = _ iz (-iJL) J e ^QW-- M ?+™)+W i + vt+Q i -yD d a i . (117) 



2Q\Q l _ -Mf + isi 

Then we substitute (117) into (116). Due to the exponentials that cut off q^, q k __ and a 1 
the integral over these variables is absolutely convergent. Therefore, one can interchange the 
integrations over 5* , q k _ and a 1 . As a result, we obtain the equality 

00 

F = lim lim lim I 0(^,^^,(3) e~^ a % (118) 

n 

where 

0(a t ,p% 1 ,(3) = (-z) n f(-z^-] x 



x / [jdYe^'[ !a ' (20 + Q '" M ' )+!(0 + J ' +<3 ' J ' h7Q + 2 ^ Q ' 2 ] . (119) 

J k Vi=0 

For the Lorentz-covariant calculation in the a-representation satisfying conditions (112), there 
is a known expression [31] 

00 

F cov = lim / y?cov(«^) e~^ a % (120) 

ae— >0 J 

n 

where 

VcoA^P 3 ) = {-i) n f ( \ x 



x 

7 



lim f TTdV e T.X i ^QXQ-- M r)+ i (QXyt+Q-vi)-^Q\ 2 -PQ 1 - 2 ] (\2i) 
,13^0 J 11 y w =o" K > 
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In Appendix 2, it is shown that in (118) the limits in 7 and (5 can be interchanged, in turn, with 
the integration over {c^}, and then with / (^—i-J^j. After that, a comparison of relations (118), 

(119) and (120), (121), clearly shows that F = F cov . Considering (116) and using Statement 1 
from Appendix 1, we obtain the equality 

F cov = lim [Udqt [ T\dq k , -Jj^L (122) 

Expression (122) differs from (114) only by the range of the integration over q k _. 

4.2 Reduction of the difference between the Light Front and 

Lorentz-covariant Feynman integrals to a sum of configura- 
tions 



Let us introduce a partition for each line, 



J dQ^ + J dQ_ = J dQ_ + (-1) J dQ 



(123) 



We call a line with integration w.r.t. the momentum Q % _ in the range (—e,e) (before removing 
the (^-functions) a type-1 line, a line with integration in the range (—00, —e) U (e, 00) a type-2 
line, and a line with integration over the whole range (—00, 00) a full line. In the diagrams, 
they are denoted as shown in Figs, la, b, and c, respectively. 



-G H- ¥r 



abode 

Figure 1. Notation for different types of lines in the diagrams: "a" is a type-1 line, "b" is a 
type-2 line, "c" is a full line, "d" is an e-line, and "e" is a H-line. 

Let us substitute partition (123) into expression (114) for F\f and open the brackets. Among 
the resulting terms, there is F cov (expression (122)). We call the remaining terms "diagram 
configurations" and denote them by Fj. Then we arrive at the relation — F cov = where 



vinB L 



and Vj is the region corresponding to the arrangement of full lines and type-1 lines in the given 
configuration. 

Note that before taking the limit in e, equations (122) and (124) can be used successfully: 
first, they are applied to a subdiagram and, then, are substituted into the formula for the entire 
diagram. This is admissible because, after the deformation of the contours described in the 
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proof of Statement 1 from Appendix 1, the integral over the loop momenta of the subdi- 
agram converges (after integration over the variables {q 1 !} of this subdiagram) absolutely and 
uniformly with respect to the remaining loop momenta {q 1 !}- Therefore, one can interchange 
the integrals over {<?+} and {q^ }. 

Thus, the difference between the LF and Lorentz-covariant calculations of the diagram is 
given by the sum of all of its configurations. A configuration of a diagram is the same diagram, 
but where each line is labeled as a full or type-1 line, provided that at least one type-1 line 
exists. 

4.3 Behavior of the configuration as e — > 

We assume that all external momenta p s are fixed for the diagram in question and 

P- t^O, EP-VO, (125) 

s' 

where the summation is taken over any subset of external momenta; all of these momenta are 
assumed to be directed inward. 

Let us consider an arbitrary configuration. We apply the term "e-line" to all type-1 lines 
and those full lines for which integration over Q- actually does not expand outside the domain 
(—re, re), where r is a finite number (below, we explain when these lines appear). The remaining 
full lines are called Il-lines. In the diagrams, the e-lines and IT-lines are denoted as shown in 
Figs. Id and e, respectively. Note that the diagram can be drawn with lines "a" and "c" 
from Fig. 1 (this defines the configuration unambiguously), or with lines "d" and "e" (then the 
configuration is not uniquely defined). 

If among the lines arriving at the vertex only one is full and the others are type-1 lines, 
this full line is an e-line by virtue of the momentum conservation at the vertex. The remaining 
full lines form a subdiagram (probably unconnected). By virtue of conditions (125), there is 
a connected part to which all of the external lines are attached. All of the external lines of 
the remaining connected parts are e-lines. Consequently, using Statement 1 from Appendix 1, 
we can see that integration over the internal momenta of these connected parts can be carried 
out in a domain of order e in size, i.e., all of their internal lines are e-lines. Thus, an arbitrary 
configuration can be drawn as in Fig. 2 and integral (124), with the corresponding integration 
domain, is associated with it. 




Figure 2. Form of an arbitrary configuration: IT is the connected subdiagram consisting of 
Il-lines, £ is the subdiagram consisting of e-lines and, probably, containing no vertices. 
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Let us investigate the behavior of the configuration as e — > 0. From here on, it is convenient 
to represent the propagator as 

y-iffw where m = 2W) or = 2Q- + to/O/ (126) 

rather than as (110). Then, in (109), M 2 = m 2 + Q^_ 2 7^ and the function / is no longer 
a polynomial. If the numerator of the integrand consists of several terms, we consider each 
term separately (except when the terms arise from expressing the propagator momentum Q l _ 
in terms of loop and external momenta). 

We denote the loop momenta of subdiagram II in Fig. 2 by q l and the others by k m . We 
make following change of integration variables in (124): 

k™ -> e k™. (127) 

Then, the integration over k™ goes within finite limits independent of e. We denote the power 
of e in the common factor by r (it stems from the volume elements and the numerators when 
the transformation (127) is made). The contribution to r from the expression 1/(2(5- +ise/Q+) 
(equation (126)), which is related to the £-line, is equal to -1. We divide the domain of inte- 
gration over k™ and q l + into sectors such that the momenta of all full lines Q\ have the same 
sign within one sector. 

In Statement 1 of Appendix 1, it is shown that for each sector, the contours of integration 
over q l _ and k™ can be bent in such a way that absolute convergence in q l + , k™, q l _ and k m 
takes place. Since, in this case, the momenta Q l _ of IT-lines are separated from zero by an 
^-independent constant, the corresponding IT-line- related propagators and factors from the 
vertices can be expanded in a series in e. This expansion commutes with integration. 

It is also clear that the denominators of the propagators allow the following estimates under 
an infinite increase in \Q+\: 



1 

2Q+Q_ — M 2 + iee 

Here c and c are e-independent constants. Note that for fixed finite Q+, the estimated ex- 
pressions are bounded as e — > 0. After transformation (127) and release of the factor ^ (in 
accordance with what was said about the contribution to r), the e-line- related expression from 
(126) becomes 

1 

2Q_ + ise/g 

where a Q+-independent quantity was used for the estimate (this quantity is meaningful and 
does not depend on e because the value of Q_ is separated from zero by an e-independent 
constant). 

We integrate first over q l + , k™ within one sector and then over q l _, k™ (the latter integral 
converges uniformly in e). Let us examine the convergence of the integral over q l + , k™ with 
canceled denominators of the e-lines (which is equivalent to estimating the expressions (129) 



< 



c\Q+\ 
1 

c e \Q A 



for LT-lines, 
for e-lines, 



(128) 
(129) 



2Q_ + iee/(Q+e) 



< 



2|Q-|' 



(130) 
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by a constant). If it converges, then the initial integral is obviously independent of e and the 
contribution from this sector to the configuration is proportional to e T . 

Let us show that if it diverges with a degree of divergence a, the contribution to the initial 
integral is proportional to e T ~ a up to logarithmic corrections. To this end, we divide the domain 
of integration over q l + , k™ into two regions: Ui, which lies inside a sphere of radius A/e (A is 
fixed), and U2, which lies outside this sphere (recall that in our reasoning, we deal with each 

sector separately). Now we estimate (128) (like (129)) in terms of - (which is admissible) 

ce\Q + \ 

and change the integration variables as follows: 

q l + ^- £ q l + , (131) 

After e is factored out of the numerator and the volume element, the integrand becomes inde- 
pendent of e. Thus, the integral converges. 

One can choose such A (independent of e) that the contribution from the domain U 2 is 
smaller in absolute value than the contribution from the domain U\. Consequently, the whole 
integral can he estimated via the integral over the finite domain U\. Now we make an inverse 
replacement in (131) and estimate (129) by a constant (as above). Since the size of the in- 
tegration domain is A/e and the degree of divergence is a, the integral behaves as e~ a (up 
to logarithmic corrections), q.e.d. This reasoning is valid for each sector and, thus, for the 
configuration as a whole. Obviously, 

a = maxav, (132) 

where a r is the subdiagram divergence index and the maximum is taken over all subdiagrams 
D r (including unconnected subdiagrams for which a r is the sum of the divergence indices of 
their connected parts). In the case under consideration, a r = uo r + + u r , where v T is the number 
of internal e-lines in the subdiagram D r . The quantities uo± are the UV divergence indices of 
the subdiagram D r w.r.t. Q±. 

Above, we introduced a quantity r, which is equal to the power of e that stems from the 
numerators and volume elements of the entire configuration. We can write r = u r _ — fi r + v r +i] r , 
where if is the index of the UV divergence in Q- of a smaller subdiagram (probably, a tree 
subdiagram or a nonconnected one) consisting of Il-lines entering D r . The term rf is the power 
of e in the common factor, which, during transformation (127), stems from the volume elements 
and numerators of the lines that did not enter D r . (It is implied that the integration momenta 
are chosen in the same way as when calculating the divergence indices of D r .) Then, up to 
logarithmic corrections, we have 

Fj ~ e\ a = mm(r, uj r _ - u T + - ^' + rf). (133) 

Consequently, for e — > 0, the configuration is equal to zero if a > 0. Relation (133) allows all 
essential configurations to be distinguished. 

4.4 Correction procedure and analysis of counter-terms 

We want to build a corrected LF Hamiltonian Hf° r with the cutoff \Q\_\ > e, which would gen- 
erate Green's functions that coincide in the limit e — > with covariant Green's functions within 
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the perturbation theory. We begin with a usual canonical Hamiltonian in the LF coordinates 
if if with the cutoff \Q % _\ > s. We imply that the integrands of the Feynman diagrams de- 
rived from this LF Hamiltonian coincide with the covariant integrands after some resummation 
[25-27]. However, a difference may arise due to the various methods of doing the integration, 
e.g., due to different auxiliary regularizations. As shown in Sec. 4.2, this difference (in the limit 
e — > 0) is equal to the sum of all properly arranged configurations of the diagram. One should 
add such correcting counter-terms to if if, which generates additional "counter-term" diagrams, 
that reproduce nonzero (after taking limit w.r.t. e) configurations of all of the diagrams. Were 
we able to do this, we would obtain the desired H™ 1 . In fact, we can only show how to seek 
the Htf T that generates the Green's functions coinciding with the covariant ones everywhere 
except the null set in the external momentum space (defined by condition (125)). However, this 
restriction is not essential because this possible difference does not affect the physical results. 

Our correction procedure is similar to the renormalization procedure. We assume that the 
perturbation theory parameter is the number of loops. We carry out the correction by steps: 
first, we find the counterterms to the Hamiltonian that generate all nonzero configurations of 
the diagrams up to the given order and, then, pass to the next order. We take into account that 
this step involves the counter-term diagrams that arose from the counter-terms added to the 
Hamiltonian for lower orders. Thus, at each step, we introduce new correcting counter-terms 
that generate the difference remaining in this order. Let us show how to successfully look for 
the correcting counter-terms. 

We call a configuration nonzero if it does not vanish as e — > 0. We call a nonzero configu- 
ration "primary" if n is a tree subdiagram in it (see Fig. 2). Note that for this configuration, 
breaking any n-line results in a violation of conditions (125); then, the resulting diagram is not 
a configuration. We say that the configuration is changed if all of the n-lines in the related 
integral (124) are expanded in series in e (see the reasoning above equation (128) in Sec. 4.3) 
and only those terms that do not vanish in the limit e — > after the integration are retained. As 
mentioned above, developing this series and integration are interchangeable operations. Thus, 
in the limit e — * 0, the changed and unchanged configurations coincide. Therefore, we always 
require that the Hamiltonian counter-terms generate changed configurations, as this simplifies 
the form of the counter-terms. Using additional terms in the Hamiltonian, one can generate 
only counter-term diagrams, which are equal to zero for external momenta meeting the con- 
dition \p s _ | < e, because with the cutoff used, the external lines of the diagrams do not carry 
momenta with \p s _ | < e. We bear this in mind in what follows. 

We seek counter-terms by the induction method. It is clear that, in the first order in 
the number of loops, all nonzero configurations are primary. We add the counter-terms that 
generate them to the Hamiltonian. Now, we examine an arbitrary order of perturbation theory. 
We assume that in lower orders, all nonzero configurations that can be derived from the counter- 
terms, accounting for the above comment, have already been generated by the Hamiltonian. 

Let us proceed to the order in question. First, we examine nonzero configurations with only 
one loop momentum k and a number of momenta q (see the notation above equation (127)). 
We break the configuration lines one by one without touching the other lines (so that the ends 
of the broken lines become external lines). The line break may result in a structure that is 
not a configuration (if conditions (125) are violated); a line break may also result in a zero 
configuration or in a nonzero configuration. If the first case is realized for each broken line, 
then the initial configuration is primary and it must be generated by the counter-terms of the 
Hamiltonian in the order under consideration. If breaking of each line results in either the first 
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or the second case, we call the initial configuration real and it must be also generated in this 
order. 

Assume that breaking a line results in the third case. This means that the resulting con- 
figuration stems from counter-terms in the lower orders. Then, after restoration of the broken 
line (i.e., after the appropriate integration), it turns out that the counter-terms of the lower 
orders have generated the initial configuration (we take into account the comment on successive 
application of equation (124); see the end of Sec. 4.2) with the following distinctions: (i) the 
broken line (and, probably, some others, if a nonsimply connected diagram arises after breaking 
the line) is not a Il-line but a type-2 line, due to the conditions \p s _\ > e; (ii) if, after restoration 
of the broken line, the behavior at small e becomes worse (i.e., a decreased), then fewer terms 
than are necessary for the initial configuration were considered in the above-mentioned series 
in e. We expand these arising type-2 lines by formula (123) and obtain a term where all of 
these lines are replaced by Il-lines or other terms where some (or all) of these lines have become 
type-1 lines. In the latter case, one of the momenta q becomes the momentum k. We call these 
terms "repeated parts of the configuration" and analyze them together with the configurations 
that have two momenta k. In the former case, we obtain the initial configuration up to dis- 
tinction (ii). We add a counter-term to the Hamiltonian that compensates this distinction (the 
counter-term diagrams generated by it are called the compensating diagrams). 

If there is only one line for which the third case is realized, it turns out that, in the given 
order, it is not necessary to generate the initial configuration by the counter-terms, except for 
the compensating addition and the repeated part that is considered at the next step. If there are 
several lines for which the third case is realized, the initial configuration is generated in lower 
orders more than once. For compensation, it should be generated (with the corresponding 
numerical coefficient and the opposite sign) by the Hamiltonian counter-terms in the given 
order. We call this configuration a secondary one. Next, we proceed to examine configurations 
with two momenta k and so on up to configurations with all momenta k, which are primary 
configurations. 

Thus, the configurations to be generated by the Hamiltonian counter-terms can be primary 
(not only the initial primary configurations but also the repeated parts analogous to them, 
called primary-like), real, compensating, and secondary. If the theory does not produce either 
the loop consisting only of lines with Q + in the numerator (accounting for contributions from 
the vertices) or a line with Q + n in the numerator for n > 1, then real configurations are 
absent because a line without Q + in the numerator can always be broken without increasing 
a (see equation (133)). It is not difficult to demonstrate that if each appearing primary, 
real, and compensating configuration has only two external line, then there are no secondary 
configurations at all. 

The dependence of the primary configuration on external momenta becomes trivial if its 
degree of divergence a is positive, the maximum in formula (132) is reached on the diagram 
itself, and a = 0. Then, only the first term is taken into account in the above-mentioned series. 
Thus, not all of the n-line-related propagators and vertex factors depend on k™ and they can 
be pulled out of the sign of the integral w.r.t. {k m } in (124). We then obtain 
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x /g rf *- n t (2Q^-^ +te) - < 134 > 

where is a domain of order e in size. Let us carry out transformations (127) and (131). For 
the denominator of the H-line, we obtain 



2(^EA; + + Ep + )(EP-)-M 2 + ^ae 2(££; + )(£p_)" v ; 

Here we neglect terms of order e in the denominator because the singularity at k™ = is 
integrable under the given conditions for a and everything can be calculated in zero order in 
e at o = 0. Thus, the dependence on external momenta can be completely collected into an 
easily obtained common factor. 



4.5 Application to the Yukawa model 

The Yukawa model involves diagrams that do not satisfy condition (112). These are displayed 
in Figs. 3a and b. We have = for diagram "a" and uj + = for diagram "b". 




Figure 3. Yukawa model diagrams that do not meet condition (112). 

Nevertheless, these diagrams can be easily included in the general scheme of reasoning. 
To this end, one should subtract the divergent part, independent of external momenta, in 
the integrand of the logarithmically divergent (in two-dimensional space, with fixed internal 
transverse momenta) diagram "a". We obtain an expression with u n < (i.e., which converges 
in two-dimensional space) and u + = 0, as in diagram "b". This means that the integral over 
q + converges only in the sense of the principal value (and it is this value of the integral that 
should be taken in the LF coordinates to ensure agreement with the stationary noncovariant 
perturbation theory). This value can be obtained by distinguishing the g+-even part of the 
integrand. 

Two approaches are possible. One is to introduce an appropriate regularization in transverse 
momenta and to imply integration over them; then, it is convenient to distinguish the part that 
is even in four-dimensional momenta q. The other is to keep all transverse momenta fixed; 
then, the part that is even in longitudinal momenta q n can be released. For the Yukawa 
theory, we use the first approach. For the transverse regularization, we use a "smearing" of 
vertices, which is equivalent to dividing each propagator by 1 + Q\ /Aj_ 2 . In four-dimensional 
space, diagram "a" diverges quadratically. Under introduction and subsequent removal of 
the transverse regularization, the divergent part, which was previously subtracted from this 
diagram, acquires the form C\ + C 2 p\. 
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After separating the even part of the regularized expression, we fix all of the transverse 
momenta again. Then it turns out that diagrams "a" and "b" in Fig. 3 meet conditions (112) 
and one can show that after all of the operations mentioned, the exponent a (see (133)) does 
not decrease for any of their configurations. Hence, they can be included in the general scheme 
without any additional corrections. 

Let us first analyze the primary configurations (see the definition in Sec. 4.4). In the numer- 
ators, appears only in the zero or one power and there are no loops where the numerators of 
all of the lines contain Consequently, one always has r > 0, fi r < 0, and rf > (see the def- 
initions in Sec. 4.3). Analyzing the properties of the expression u r _ — uo r + for the Yukawa model 
diagrams, we conclude from (133) that a > always holds. The general form of the nonzero 
primary configurations with a = is depicted in Fig. 4. Note that they are all configurations 
with two external line. 
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Figure 4. Nonzero configurations in the Yukawa model: p is the external momentum, and 7 + 
or 7~ symbols on the line indicate that the corresponding term is taken in the numerator of 
the propagator. In configuration "b", the part that is proportional to 7 + is taken. 

Further, it is clear that there are no nonzero real configurations (see the comment at the 
end of Sec. 4.4), and it can be shown by induction that there are no nonzero compensating or 
secondary configurations either (the definitions are given in Sec. 4.4 also). Thus, only primary 
or primary-like configurations can be nonzero and all of them have the form shown in Fig. 4. It 
can be shown that their degree of divergence a is positive and the maximum in formula (132) 
is reached for the diagram itself. Thus, the reasoning above and below formula (134) applies 
to them. Then, denoting the configurations displayed in Figs. 4a-d by D a - D d , we arrive at 
the equalities D a = ^-C a , D b = ^-C b , D c = C c and D d = C d , where the expressions C a - C d 
depend only on the masses and transverse momenta, but not on the external longitudinal 
momenta, and have a finite limit as e — > 0. 

Now we assume that D a - D d are not single configurations but are the sums of all configu- 
rations of the same form and that integration over the internal transverse momenta has already 
been carried out, (with the above-described regularization) . In four- dimensional space, the di- 
agrams D a and D b diverge linearly while D c and D d diverge quadratically. Therefore, because 
of the transverse regularization, the coefficients C c and C d in the limit of removing this regu- 
larization take the form C\ + C 2 p±, where C\ and C2 do not depend on the external momenta 
(neither do C a , C b )). Thus, to generate all nonzero configurations by the LF Hamiltonian, only 
the expression 

H c = C lV 2 + C 2 p\ v? 2 + C 3 ^ ^ V, (136) 
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should be added, where f and ip are the boson and fermion fields, respectively, and Q, are the 
constant coefficients. 

Comparing (136) with the initial canonical LF Hamiltonian, one can easily see that the 
found counter-terms are reduced to a renormalization of various terms of the Hamiltonian (in 
particular the boson mass squared and the fermion mass squared without changing the term 
linear in fermion mass). The explicit Lorentz invariance is absent but results are Lorentz 
invariant because the counterterms compensate the violation of Lorentz invariance inherent to 
chosen LF formalism. 

Note that the second approach, mentioned at the beginning of this section, can give the 
same results. The only difference is that in two-dimensional space, the contributions from the 
configurations displayed in Fig. 3 would additionally depend on external transverse momenta. 
However, this dependence disappears after integration over internal transverse momenta with 
the introduction and subsequent removal of an appropriate regularization. 

In the Pauli Villars regularization, it is easy to verify that the expression u r _ — uo r + — fi r + rf 
from (133) increases. This is because the number of terms in the numerators of the propagator 
increases. Then, the contribution from the e-lines does not change, while the n-lines belonging 
to D r make zero contribution to uo r _ — uo r + and rf , but —1 contribution to fi r . Since r > 0, this 
regularization makes it possible to meet the condition o > for the configurations that were 
nonzero (one additional boson field and one additional fermion field are enough). Then it turns 
out that the canonical LF Hamiltonian need not be corrected at all. 

Obtained results agree with the conclusions of the work [25], where a comparison of LF 
and Lorentz-covariant methods was made for self-energy diagrams in all orders of perturbation 
theory and for other diagrams in lowest orders. 



4.6 Application to QCD(3+1) 

Applying the LF Hamiltonian approach to gauge theories under regularization (1), where zero 
modes of fields are thrown out, one has to use the gauge = (see, for example, sect. 2.3). To 
carry out successfully renormalization procedure for this scheme it is necessary to take gauge 
field propagator in accordance with Mandelstam-Leibbrandt prescription [28, 29] (which allows 
to perform Euclidean continuation, see [32-34]). Such a propagator has the form 

( gm _ ( m^2 Q+ ) = L, - q»"»,tQ»>. 2 iQ n- } ) , (137) 



Q 2 + iaeV 2Q+Q_+iae ^ + J Q 2 + iae \ 2(Qn*)(Qn) + 

where n + — 1, n_ = n± — 0, n*_ = 1, n* + = n*^ = 0. 

The formalism described in sect. 4.1-4.4 was such that it could be applied to a theory 
with the propagator (137) (at fixed transverse momenta Q± ^ 0). It turns out that there are 
nonzero configurations with arbitrarily large numbers of external lines. An example of such a 
configuration is given in Fig. 5. 

Indeed, using formula (133), we can see that for the configuration in Fig. 5, r = and, thus, 
er < 0, i.e., this is a nonzero configuration. It is also clear that introduction of the Pauli- Villars 
regularization does not improve the situation because it does not affect r. 

The difficulty is that the distortion of the pole in (137) due to LF cutoff \Q-\ > e > 
does not disappear in the limit e — > 0, and infinite number of new counterterms are required 
to compensate this distortion [11]. The simplest way to avoid this difficulty is to add small 
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Figure 5. Nonzero configuration with an arbitrarily large number of external lines in a gauge 
theory. The symbols 7- 1 - on the lines and the symbols + or _L the vertices indicate that the 
corresponding terms 7 + or 7- 1 - are taken in the numerators of propagators and in the vertex 
factors. 

mass-like parameter fj 2 in the denominator: 



2Q + Q_+ias 2Q + Q_ - /i 2 + iae' 

and take the limit e — > before // — > 0). To describe this modification with local Lagrangian 
we need to introduce ghost fields A' in addition to conventional A^. We write the free part 
of pure gluon Lagrangian as follows (using higher derivatives and the parameter A for UV 
regularization): 

a> = -\[r^ (1 + 5) - (1 + % + ^=)fs) > ( 139 ) 

where /£, = d^A a v - d u A®, = d^AJ - d v A'£ and A a _ = A m = 0. Interaction terms depend 
only on summary field A a ^ = A a ^ + A'^. 

At fixed ji and A we get a theory with broken gauge invariance but with preserved global 
S'£/(3)-invariance. We put into the Lagrangian all necessary interaction terms (but with un- 
known coefficients) including those that are needed for UV renormalization: 

L = L + c^A^A^ + c m d ll A a Xn* a d a A a ^ + c^A^ d v A a,v + 

+c 11 n^n* a d fl A a ' a d u A a ' u + c l2 n^n* a d II A a ' a n u n^d u A^ + c 2 A^A a '» + 
+c 3 f abc A^A b iy d^A c '' / + c 31 f abc A a ^Ain a n*>*d a A c ' u + 
+A*A b u A;A d s (cJ abe f cde g^g' /S + 5 ab 5 cd (c^g uS + c 6 ^V 5 )) . (140) 

These terms are local and can be taken in Lorentz covariant form due to the restoring of this 
symmetry in the e — > 0, \x — > limit. However in constructing of these terms both vectors n M 
and n* u , included in the definition of the propagator (137), can participate. 

For this theory one can apply the formalism described in sect. 4.1-4.4 to compare (at e — > 0) 
the LF perturbation theory and that taken in Lorentz coordinates within the same regulariza- 
tion scheme. We find that the difference between mentioned perturbation theories can be com- 
pensated by changing of the value of coefficient c 2 before the term of gluon mass form A a ^A^> a in 
naive LF Hamiltonian of this theory. After that we can analyse further our regularized theory 
in Lorentz coordinates and even make Euclidean continuation. 
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It is possible to prove by induction to all orders of perturbation theory that in the limit 
fi — > 0, A — > oo our theory can be made finite and coinciding with the usual renormalized 
(dimensionally regularized) theory in Light Cone gauge [32, 33] (for all Green functions). This 
was done in the paper [14], but there the terms of the Lagrangian, including the vector rf 1 , 
were missed. Right expression for LF Hamiltonian should correctly take into account the 
contribution of all terms, written in the (140). 

The values of the unknown coefficients q before all counterterms in (140) must be chosen 
so that the Green functions in each order coincided (after removing the regularization) with 
those obtained in conventional dimensionally regularized formulation and therefore satisfied 
Ward identities. Besides, we need to correlate the limits /i — > and A — > oo to avoid infrared 
divergencies at /i — > 0. It is sufficient to take /i = /i(A) and to require that /iA — > and 
Oog'/O/A - 0. 

Our resulting LF Hamiltonian for pure SU(3) gluon fields contains 11 unknown coefficients, 
including coefficient before gluon mass term that takes also into account the difference between 
LF and Lorentz coordinate formulations of our regularized theory. The generalization of our 
scheme for full QCD with fermions is described in [14]. In this case there are 20 unknown 
coefficients in the LF Hamiltonian. We hope that it is possible to find an analog of Ward 
identities relating the coefficients q at fixed A. This problem seems very important for our 
approach. 



4.7 Application to QED(1+1) 

The procedure, described in previous section, allows to construct LF Hamiltonian for QCD 
in four-dimensional space-time, but such LF Hamiltonian contains many additional fields and 
unknown coefficients. This complicates calculations with this Hamiltonian. Moreover, because 
only the perturbation theory with respect to the coupling constant was analyzed, there could 
remain purely nonperturbative effects that are not taken into consideration. It is therefore 
useful to consider an example of "nonperturbative" (with respect to usual coupling constant) 
construction of LF Hamiltonian for gauge theory, which is possible for two-dimensional quantum 
electrodynamics (QED(1 + 1)), i.e. for "massive" Schwinger model. 
The QED(1+1), defined originally by the Lagrangian 

L = —F^F^ + ^(i! m D^ - M)tt, (141) 
can be transformed to its bosonized form [35, 36], described by scalar field Lagrangian 

L = ± (d„<S>d»<S> - m 2 $ 2 ) + — cos(# + v 7 !^) , (142) 

where m = e/^ is the Schwinger boson mass (the e is original coupling), C = 0.577216 . . . 
is the Euler constant, and the is the "^"-vacuum parameter, which takes into account the 
nontriviality of QED(1+1) quantum vacuum due to instantons. Here the fermion mass M 
plays the role of the coupling in bosonized theory so that perturbation theory in this coupling 
corresponds to chiral perturbation theory in QED(1 + 1). The nonpolynomial form of scalar 
field interaction leads in perturbation theory to infinite sums of diagrams in each finite order. 
It can be proved [16, 37] that some partial sums of these infinite sums are UV divergent in 
the 2nd order, whereas for full (Lorentz-covariant) Green functions these divergencies cancel 
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(remaining only the divergent vacuum diagrams). Therefore physical quantities are UV finite 
in this theory. Only at intermediate steps of our analysis we need some UV regularization. 

We compare LF and Lorentz-covariant perturbation theories for such bosonized model us- 
ing an effective resummation of perturbation series in coordinate representation for Feynman 
diagrams [16, 17] and also using the formalism described in sect. 4.1-4.4. The results of this 
comparison can be formulated as follows. 

The difference between considered perturbation theories can be eliminated in the limit of 
removing regularizations if we use instead of the naive LF Hamiltonian 

H = fdx~ (^m 2 : y? 2 : - 1 - e w : e* : e' ie : ■) , 
J V87T 2 2/ 

Mme c /— — . 

7 = — , y? = V47r$, \p-\>e>0, (143) 

2n 

the "corrected" LF Hamiltonian: 

H = Jdx- (— m 2 : V 2 :-B: e itp : -B* : e~ i(p :) - 

_ 27re - 2C ^ / dx'dy- (: e ^(*-) e -*^-) : -l) e(\x~ - y~\ - a) V ^ X ~ ~ ■ (144) 
7 v 7 |x~ — y~| 

Here the terms, linear in B and 5* (new coupling constants), are of the same form as in naive 
Hamiltonian; only the term, containing the \B\ 2 , is of new form (nonlocal in x~). The a is 
the UV regularization parameter, and the v(z) is some arbitrary continuous function rapidly 
decreasing at the infinity and going to unity as z — > 0. The coupling B can be perturbatively 
written as a series in 7: 

B = ^ + E 7 n ^n- (145) 

2 n=2 

On the other side, it is related to the sum of all connected "generalized tadpole" diagrams 
(i.e. diagrams with external lines attached to only one vertex), which is described by the 
"condensate" parameter A = ^(f2| : e l ^ +9 ^ : of the Lorentz-covariant formulation (the 
is the physical vacuum state in this formulation): 

B + \B\ 2 w = A, (146) 



27re _2C r , 9(\x - ea) . s 

w = r— / rfa;~ , , -ufx ). (147) 

rrr 7 |rr _ | 

The equation (146) can be solved with respect to the B: 

B = ~ + J-^- + — - A" 2 + iA", (148) 
2w V 4w 2 w 

where A = A' + iA", and the sign before the root respects the perturbation theory. Within 
the perturbation theory in 7 one cannot remove UV regularization (i. e. to put a — > and 
therefore w — > 00) in this expression due to UV divergencies of the coefficients B n . However, 
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taking into account the validity of the equation (146) to all orders in 7, we can consider it 
beyond the perturbation theory. Then we use the estimation for the A at a — > [16]: 

7 2 

A = —w + const (149) 
4 

and get for the B in a — > limit UV finite result: 



£ = sign(cosfl)y ^ - A" 2 + L4" = |e^, (150) 

so that all information about the condensate is contained in the phase factor e ld : 

Tm A 

sin 9 = 2 = (n\ : sin(^ + 9) : (151) 

7 

Then we can make a transformation, inverse to the bosonization, but on the LF. Actually 
we need the expression only for one independent component ip + of the bispinor field (^ + J 
due to the LF constraint, permitting to write the in terms of ip + . One can use the exact 
expression for the ip + in terms of the tp obtained in the theory on the interval \x~\ < L with 
periodic boundary conditions [36, 16]. We need only to modify our corrected bosonized theory 
by using discretized LF momentum p_ instead of continuous one and hence replacing the cutoff 
parameter e by tt/L. The necessary formulae for the tp + has the following form [16, 36] (we 
choose here antiperiodic boundary conditions for the fermion fields): 

rl)+{x) = -Le^e^^^e^ 1 " : e'^ : . (152) 
v 2L 



The operator u and the charge operator Q are canonically conjugated so that the ip + defined 
by the equation (152) has proper commutation relation with the charge. On the other side the 
operator e luJ shifts Fourier modes ip n of the field ip + [16, 36]: 

e^e-*" = ^ n+ i. (153) 

If we separate the modes related with creation and annihilation operators on the LF putting 

iM*) = b n e~^-k)*- + £ d+e^^-Y 6 n |0) = d n |0) = 0, (154) 

V 2L \ n > x „> / 

we can define the operator e tuJ uniquely by specifying its action on the LF vacuum |0) as follows: 

e-|0) = fc+|0), e — 10) = 4|0). (155) 

In such sense this operator is similar to a fermion. 

We can now rewrite our corrected boson LF Hamiltonian in terms of tp + and e lw . The result 
is remarkably simple: 

H= /^-(^(9: 1 [^ + ]) 2 -^9:V + - (^e-*e-d+ + h.c)). (156) 
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This fermionic LF Hamiltonian differs from canonical one (in corresponding DLCQ scheme) 
only by last term, depending on zero modes and vacuum condensate parameter 9 which can be 
related to chiral condensate by transforming the variables in the equation (151): 

2vr 3 / 2 

sin# = : ^ 7 5 ^ : \Q). (157) 

e e° 

Let us remark that the presents of linear in M term in LF Hamiltonian (156) can be 
considered as a results of a modification of the LF constraint, connecting the ip- with the 
An analogous modification of this constraint was got in the paper [38] where the method of 
exact operator solution of massless Schwinger model was applied. 

The constructed LF Hamiltonian (156) was applied to nonperturbative numerical calculation 
of mass spectrum of QED(1+1) [18]. The results of this calculation were compared with those 
of lattice calculations in Lorentz coordinates [39, 40]. 

The calculations were carried out in wide domain of values of fermion mass M for all values 
of the parameter 9, which is a function of the M/e and vacuum parameter 9. We do not know 
exactly this function, but know that it must be zero at 9 = and be equal to n at 9 = ir. 

For 9 = the obtained spectrum is bounded from below at any values M, and is in good 
agreement with the results of the paper [39] for two bound states of lowest mass. 

For the value 9 = tt, at which phase transition can take place, the obtained spectrum for 
the lowest bound state agrees well with the results of the paper [40] for sufficiently small M; 
at greater M we start to see the disagreement with that paper, and then at M, greater some 
critical value, the spectrum becomes unbounded from below. This critical value approximately 
coincides with the point of phase transition found in [40]. It can be supposed that at M greater 
than the point of phase transition the calculations with our LF Hamiltonian, constructed via 
the analysis of perturbation theory in M, become incorrect. 

Our calculations show that LF Hamiltonian (156) can give good results only in the limited 
domain of the parameters: from perturbative domain to some values which define the limits of 
applicability of our LF Hamiltonian. In the domain where our Hamiltonian becomes incorrect 
we need to take into account nonperturbative (in M/e) effects. Such an investigation can be 
useful for finding a LF Hamiltonian approach to realistic gauge theories. 

The LF Hamiltonian (156) includes the operator ou, which has no simple expression in terms 
of field operators. It is defined only by its properties (153), (155). Due to this fact the expression 
for the Hamiltonian depends essentially on the form of the regularization, i.e. \x~\ < L and 
antiperiodic boundary conditions in x~ for the field ip. Now we have found a possibility to 
rewrite the expression (156) in such a way that it contains only fermion field operators, and 
describes in the limit of removing the regularizations the same theory as the Hamiltonian (156). 
This new expression has at 9 = 9 = the following form: 

H=f dx- (cT 1 ^]) 2 + (4* + bfh) - ^VcTV) • (158) 

Preliminary calculations of the mass spectrum, produced by this Hamiltonian, show that in the 
limit of removing the regularization, L — > oo, results indeed coinside, with a good accuracy, 
with those for the bound state mass spectrum found here for the Hamiltonian (156). Work 
on studing of the Hamiltonian (158) spectrum and also on the constructing the analogous 
expression for the Hamiltonian at 9 ^ will be continued in future. 
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5 Transverse lattice regularization of Gauge Theories 
on the Light Front 

The introduction of space-time lattice for gauge-invariant regularization of nonabelian gauge 
theories is well known [41]. Gauge invariant regularization in continuous space-time is also 
known [42] but it seems not suitable for the LF quantization. For the LF formulation only 
the lattice in transverse coordinates x x ,x 2 is used. In this formulation it is convenient to 
define variables so as to have the action polynomial in these variables [43, 44, 19, 20]. Such 
a regularization is not Lorentz invariant, and one can only hope that Lorentz invariance can 
be restored in continuous space limit. Nevertheless many attempts to apply LF Hamiltonian 
formulation with periodic boundary conditions in x~ , combined with transverse space lattice, 
are undertaken (for "color dielectric" type models [45-48]). In all of these works zero modes of 
fields are thrown out, so that, in fact, gauge invariance is violated. 

We consider canonical LF formulation of gauge theories, regularized in gauge-invariant 
way. To achieve this goal we introduce transverse space lattice, discretize the momentum 
p_ according to the prescription (b) (see introduction, equation (2)) with all zero modes of 
fields included and apply the so called "finite mode" ultraviolet regularization in p_. The last 
means a cutoff in eigenvalues of covariant derivative operator £L in the expansion of lattice field 
variables in eigen functions of this operator. These field variables are lattice modification of 
transverse components of usual gauge fields. They are described by complex matrices, defined 
on lattice links. Only such variables admit mentioned above "finite mode" regularization (for 
fermion fields analogous method was applied in [49, 50]). 

It is interesting that in the framework of this formulation one can avoid complicated canon- 
ical 2nd-class constraints, usually present in canonical LF formalism in continuous space. This 
greatly simplifies canonical quantization procedure. However the absence of explicit Lorentz 
invariance of the regularization scheme makes the investigation of the connection with conven- 
tional Lorentz-covariant formulation difficult. In particular, there is a problem of the description 
of quantum vacuum as common lowest eigen state of both operators P_ and P + . 

5.1 Gauge-invariant action on the transverse lattice 

At first we introduce particular ultraviolet regularization via a lattice in transverse coordinates 
x x ,x 2 and choose variables so as to have the action, which is polynomial in these variables 
[43, 44]. Furthermore, we use the described in the introduction gauge- invariant regularization 
(b) (see equation (2)) of singularities at p_ — > and gauge-invariant ultraviolet cutoff in modes 
of covariant derivative £L (then ultraviolet regularization of the theory is complete). For 
simplicity further we consider again the U (N) theory of pure gauge fields because this example 
is technically more simple than the SU (N) theory. 

The components of gauge field along continuous coordinates x + , x~ can be taken without 
a modification and related to the sites of the lattice. Transverse components are described by 
complex N x N matrices M^(x), k = 1,2. Each matrix M^(x) is related to the link directed 
from the site x — to the site x. The transverse vector connects two neighbouring sites on 
the lattice being directed along the positive axis x k (\ek\ = a), see fig. 6. The matrix M£(x) is 
related to the same link but with opposite direction, see fig. 7. In the following the usual rule 
of summation over repeated indices is not used for the index k. If this summation is necessary, 
the symbol of a sum is indicated. 
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M k {x) 



x - e k 



x 



axis x 



Figure 6. 



M+(x) 



x - e k 



x 



axis 



x 



Figure 7. 



The elements of these matrices are considered as independent variables. This makes the 
action polynomial. For any closed directed loop in the lattice we can construct the trace of the 
product of matrices M k {x) sitting on the links in the loop and order these matrices from the 
right to the left along this loop. For example the expression 



Tr {m 2 (x)M 1 (x - e 2 )M+(x - ei )M+(x)} 
is related to the loop shown in fig. 8. 



(159) 



x — ei 



Mt{x-e 1 ) 




x - ei - e 2 



M±(x-e 2 ) x-e 2 



axis x 



A 



axis x 



Figure 8. 

It should be noticed that a product of matrices related to closed loop, consisting of one and 
the same link passed in both directions, is not identically unity because the matrices M k are 
not unitary (see, for example, fig. 9). 

The unitary matrices U (x) of gauge transformations act on the M and M + in the following 

way: 



M k (x) -> M' k (x) = U(x)M k (x)U + (x - e k ), 



(160) 



M+(x) -> M'+(x) = U(x - e k )M+(x)U + (x). 



161) 
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Figure 9. 



A trace of the product of the matrices, related to closed loop along lattice links, is invariant 
with respect to these transformations. To relate the matrices M k with usual gauge fields of 
continuum theory let us write these matrices in the following form: 



M k (x) = 1 + gaB k {x) + igaA k (x), B+ = B k , A+ = A k . 



(162) 



Then in the a — > limit the fields A k (x) coincide with transverse gauge field components, 
and the B k (x) turn out to be extra (nonphysical) fields which should be switched off in the 
limit. Below we show how to get this. 

The analog of the field strength 



d^A u - dvAp - ig[An, A u ], 



multiplied by i, can be defined as follows: 

G + _ = iF + _, F + _(x) = d + A_(x) - d_A + {x) - ig[A + (x), A_(x)}, 

G± >k (x) = — [d±M k (x) - ig (A±(x)M k (x) - M k (x)A±(x - e k ))\ , 
ga 

G 12 (x) = -\ [M 1 (x)M 2 (x - e x ) - M 2 (x)M 1 (x - e 2 )\ . 

ga 

Under gauge transformation these quantities transform as follows: 

G + _(x) -> G' + _(x) = U(x)G + _(x)U + {x), 
G ± , k (x) -> G' ±jk (x) = U(x)G±, k (x)U + (x - e k ), 
G 12 (x) -> G' 12 (x) = U(x)G 12 (x)U + (x - ei - e 2 ). 

We choose a simplest form of the action having correct naive continuum limit: 

L 



T,J dx+ J dx ~ Tr G\_G + _ + £ {G + +k G_ k + G + _ k G +k ) - Gf 2 G 12 

v-L t _ k 



(163) 



(164) 



(165) 



+ S rn , (166) 



where the additional term S m gives an infinite mass to extra fields B k in the a — > limit: 



S m — 



m 2 (a) 



4gi 



J2 f dx+ [ dx'^Tr (M k + (x)M k (x) -I 



— ■> -m 2 (a) J Sx L J dx + J dx~ Tr { B l) > 



m(a) — > oo. (167) 

a— >0 



It is supposed that this leads to necessary decoupling of the fields B k . 
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5.2 Canonical quantization on the Light Front 

Let us fix the gauge as follows: 

<9_A_ = 0, A^x) = 5 ij v j (x ± ,x + ). (168) 

For simplicity below we denote the argument of quantities, not depending on the x~ , again by 
x. Let us remark that starting with arbitrary field A^, periodic in x~ , it is not possible to take 
zero modes of the A_ equal to zero without a violation of the periodicity. But it is possible to 
make the A_ diagonal as in the equation (168) [21-24]. 
Then the action (166) can be written in the form: 

S = a 2 J2j dx + J dx-\j2 [zFl_(x)d + v\x)\ + 

x-L _ L I i 

^EE + + E4WQ J, W - w(4 (169) 



(gay 



i,j k ' i,j 



where 



D_M l k 3 (x) = [d- - igv\x) + igv 3 (x - e k )) M l k J (x), 
D_MH + (x) = (&_ + igv\x) - igv j {x - e k )) M l k j+ (x), 

D_F^_{x) = (d- - igv\x) + igv j {x)) F%_(x), (170) 

the A\ (x) play the role of Lagrange multipliers, 

Q^ i (x) = 2D_Fl i _(x) + 



+^EE M? + (x)D-Ml j '(x) - M^ + {x + e k )D.M(\x + e k ) 
ga k L 



0, (171) 



- (p-M%' + M»'(x) + [D^ + (x + e fc )) M(\x + e k ) 
are gauge constraints and 

W = E (Ft + Ft + G% + G%) + H m (172) 



is Hamiltonian density. The term 7i rn can be obtained from the expression (167) in standard 
way. 

The constraints can be resolved explicitly by expressing the F+_ in terms of other variables, 
but zero mode components -^+_(o) can n °t be found from constraint equations and play the 
role of independent canonical variables. Zero modes QY Jx ± ,x + ) of the constraints remain 
unresolved and are imposed as conditions on physical states: 

Q% ) (x L ,x + )\* phy9 ) = 0. (173) 
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In order to find complete set of independent canonical variables we write Fourier transformation 
in x~ of fields M k (x) in the following form: 



M*(x) = -f= J2 {Q(pn + gv\x)-gv^x-e k ))Mii k {x\x + ) + 



n=—oo 



+© (~Pn ~ gv\x) + gv j (x - e k j) M t n j k + (x ± ,x + )} x 



X 



p n + gv\x) - gv j (x-e k ) V ' e~ iPnX ~ , (174) 



where 



„ , , f 1, p > 7T 

^ = \0, p<0 ' Pn = jn, neZ. (175) 

Due to the gauge (168) this Fourier transformation coincides with the expansion in eigen modes 
of the operator £)_ . Therefore the ultraviolet cutoff in these modes, which we will apply, reduces 
to the following condition on the number of terms in the sum (174): 

7T 

\p n + gv % (x) -gv 3 (x-e k )\ < —n, (176) 

where the n is integer parameter of ultraviolet cutoff. Let us stress that this regularization is 
gauge invariant. 

The action can be rewritten in the following form (up to nonessential surface terms): 



S = a 2 J2[dx+ |]T ALF i ;_ {0) d + v i + 



+i E E E'M^+Kl + 2LJ2 A« {0) Qf 0) - H(x) \ , (177) 

" i,j k n i 



where the J2' n means that the sum is cut off by the condition (176), and the 7i is obtained from 
the TC via the substitution of the expression 

Ft = (Ft - fi»F*_ i0) ) + S ij Ft {0) , (178) 

where the F%_ — 5 lj F™_^ are to be written in terms of the M% k , M l J k + , v l by solving the 

constraints (171) and using the equation (174). The -F+_( ) remain independent. The (?l 2 are 

also to be expressed in terms of the M% k , M^ k + , v l via the equations (164), (174). 
We have the following set of canonically conjugated pairs of independent variables: 

{v\ ff = 4La 2 Ft ( o)}, {Mt *M^ k + ) . (179) 

In quantum theory these variables become operators which satisfy usual canonical commutation 
relations: 



■'x J -,x'- L - ) 



[v\x),W(x')] x+=0 = t^5 x 
[M l n j k (x), M% + (x% +=Q = 5 ti '5"'5 nn ,5 kk ,5 x ^ x „; (180) 
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the other commutators being equal to zero. 

Let us remark that the condition (168) does not fix the gauge completely. In particular, 
discrete group of gauge transformations, depending on the x~ , of the form 

U l n j (x) = F'exp ji^Or^x- j , (181) 

where n^x -1 ) are integers, remains, and, of course, transformations, not depending on the x~ , 
are allowed. Under the transformations (181) canonical variables change as follows: 

v\x) — ► v\x) - n\x^), IT — > IT, 

Ki(^) — > K'ki^l n' = n + n\x^) - n^x 1 - e k ). (182) 

Let us write the expression for quantum operators QJ Q - ) (x ± ,x + ), which define the physical 
subspace of states. We fix the order of the operators in such a way as to relate with classical 
expression G+ l/ G lxv a quantum one of the form: 

\(0+G^ + G^G+). (183) 

We remark that other choices of the ordering do not admit reasonable vacuum solution. Then 
the operators Qf Jx ± : x + ) have the following form in terms of canonical variables: 

2LQf 0) (x ± , x+) = -j^EEE' [ £ (Pn + 9v j (x + e k ) - gv\x)) x 

40 3 k n 

x (M^ k + (x + e k )Mg(x + e k ) + M^x + e k )M^ k {x + e k )) - 

-e (p n + gv\x) - gv\x - e k )) (m^ \x)M l n j k (x) + M^x) \x))] , (184) 

where 

4P) = {X III (185, 
One can easily construct canonical operator of translations in the x~: 
P- an = \HHH E'Pn e (p n + gv\x) - gv\x - e k )) x 

x (Mg + (x)Mg(x) + Mg{x)M^ + (x)) . (186) 

This expression differs from the physical gauge invariant momentum operator P_ by a term 
proportional to the constraint. The operator P_ is 

2 L 

p - = yEE / dx~ Tr (G + _ k G- k + G. k G± k ) = P™' + ALa 2 £ £ v'Qfa = 



rA k 



= 2 E E E E' W + - 9v 3 {x - e k ) 



i,j k n 



X 



x (M^ k \x)M^{x) +Mi l l( X )MH k + (x)) = 
EEEE' \Pn + 9v\x) - gv\x - e fc )| (m^ + (x)M^(x) + ±) . (187) 
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Let us choose a representation of the state space, in which the variables v l (x) are the multi- 
plication operators. The states are described in this representation by normalizable functionals 
F[v] of classical functions v l (x) (in fact by functions, depending on the values of the v l in 
different points x x due to the discreteness of these x -1 ). One can define full space of states as 
direct product of Fock space, in which the M% k (x) and M^ k (x) play the role of creation and 
annihilation operators, and the space of functionals F[v]. Let us call M- vacuum the states of 
the form |0) • F[v], where the |0) satisfies the condition 

M% t (x)\0)=0, (188) 

and the F is some functional. Arbitrary state can be represented in the form of linear combi- 
nation of vectors \{m}\ F) of the form 

nnnn'(^ + (-)) m " fcW io)^M 

x-L i,j k n 

with different nonnegative integer functions rn^x^) and functionals F. One can define or- 
thonormalized set of such functionals if necessary. 

One can see from (187) that the state, corresponding to the absolute minimum of the P_ 
must satisfy the conditions (188), i.e. to be a M- vacuum. The value of the P_ in this state can 
be written in the form 

<0;P|P_|0;P> = 

= \ I II II EEEE' \Pn + gv'ix) - gvi(x - e k ) \ \F[v}\ 2 . (190) 

x ± i x 1 - i,j k n 

Remind that the J2' denotes the sum in n limited by the condition (176). If one uses this 

n 

condition and shifts the index n in these sums by integer part of the quantities (gL(v % (x) — 
v\x — e k ))/ir), one sees that the dependence on the v l cancels in sums over n and that the 
expression (190) does not depend on the F[v] if it is normalized to unity. Thus the momentum 
P_ has the minimum in all M-vacua. One can make the value of the P_ in these vacua equal 
to zero by subtracting corresponding constant from the operator P_. 

Let us show that M-vacua are the physical states, i. e. satisfy the condition (173). Indeed, 
in M-vacua this condition looks as follows: 

E E E' [ £ (Pn + 9v j (x + e k ) - gv\x)) -e(p n + gv\x) - gv j {x - e k j)] F[v] = (191) 

j k n 

and is satisfied for any F[v], because for every link in the sum (191) the numbers of positive 
and negative values of the ^-functions are equal. For arbitrary basic states (189) analogous 
conditions have the following form: 

EEE' [ £ (Pn + gv J (x + e k ) - gv l {x])m^ k {x + e k )- 

j k n 

- e (p n + gv\x) - gv j (x - e k )) m^ k {x)] F[v] = 0. (192) 
The eigen values p- of the operator P_ can be found from the equation 

Y,Y,Y,Y,'\Pn + gv'ix) - gv\x - e k ) | m%{x)F[v\ = p.F[v] (193) 
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where the functional F[v] is normalized. 

To define physical vacuum state correctly one must consider not only states, corresponding 
to the minimum of the operator P_, but also to the minimum of the operator P + . One can try 
to do this via minimization of the P+ on the M- vacua, i. e. on the set of states with p_ = 0. The 
expression (0|P+|0), where |0) is the Fock space vacuum w. r. t. the M n k, M n k + , depends on 
the functions v l (x) (which enter, in particular, into the "normal contractions" of the operators 
Mfc, M£) and on the operators II*, canonically conjugated to the v l (x). The expectation value 
of this expression is to be minimized on the set of functionals F[v]. The resulting functional 
F[v] must decrease in the vicinity of those values of the v l (x), for which the operator £L has 
zero eigen value, because the Hamiltonian is singular at these values (it is seen, for example, 
from the expansion (174)). 

The vacuum state constructed in such a way strongly differs from the usual vacuum in 
continuous space theory, because for M-vacuum we get zero expectation values of the operators 
Mfc, but not of the operators (Mk — I) /go, = Bk + iAk, related to usual gauge fields. Beside 
of this the condition of the unitarity of the matrices in the continuum limit (or equivalent 
condition of switching off the nonphysical fields B k ) cannot be fulfilled in such a vacuum. This 
disagreement with conventional theory is caused by the absence of explicit Lorentz invariance 
in our formulation, that leads to different quantum states, corresponding to the minima of the 
operators P_ and P + . It is not clear whether these states can be made coinciding at least in 
the limit of continuous space. This requires further investigation. 

Nevertheless the method of the quantization of gauge theories on the LF, described here, 
can be useful for completely gauge-invariant formulation of some effective models, based on 
analogous formalism (but without complete gauge invariance due to throwing out of all zero 
modes of fields and due to the absence of gauge-invariant regularization of ultraviolet divergen- 
cies in the p_). Such models are described, for example, in papers [45, 46], where the ideas of 
papers [47, 48] were developed. 

6 Conclusion 

The LF Hamiltonian approach to Quantum Field Theory briefly reviewed here is an attempt 
to apply a beautiful idea of Fock space representation for quantum field nonperturbatively in 
the framework of canonical formulation on the LF. The problem of describing the physical 
vacuum state becomes formally trivial in this approach because such vacuum state coincides 
with mathematical vacuum of LF Fock space. 

However the breakdown of Lorentz and gauge symmetries due to regularizations generates 
difficulties in proving the equivalence of LF formalism and the usual one in Lorentz coordinates. 
This problem can be solved for nongauge theories but turns out to be very difficult for gauge 
theories in special (LF) gauge which is needed here. Nevertheless we hope that these difficulties 
can be overcome by finding a modified form of canonical LF Hamiltonian which generates the 
perturbation theory equivalent to usual covariant and gauge invariant one. 

Acknowledgments. The work was supported by the Russian Foundation for Basic Re- 
search, grant no. 05-02-17477 (S.A.P. and E.V.P.). 
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Appendix 1 



Statement 1. If conditions (112) are satisfied, then, for fixed external momenta p s and p s __ ^ 
Vs 7 the equality 



limlim T\dq k + T\dq k H % ',nr Ftv IT2 \ 
^07-0 J 11 y + i 11 y - n .( 2 g^_ - + iae) 



V \4 V ruW + e - mj> + i8e ) 1 ' 

/ioWs tu/ii/e t/ie expressions appearing in (A. 1.1) exist and the integral over on the right- 
hand side is absolutely convergent. It is assumed that the momenta of lines Q l are expressed in 
terms of loop momenta q k , V e is the domain corresponding to the presence of full lines, type-1 
lines, and type-2 lines (the definitions are given following formula (123)), Bl is the sphere of 
radius L, where L > S max \p s _\, and S is a number depending on the diagram structure. 

Let us prove the statement. For each type-1 line in (A. 1.1), we perform the following 
partitioning: 



t 1 —t CXJ 

J dQl = J dQ l _ + (-1) J dQl + J dQ[ 



(A.1.2) 



Then both sides of equation (A. 1.1) become the sum of expressions of the same form in which, 
however, the domain V £ corresponds to the presence of only full and type-2 lines. It is clear that 
by proving the statement for this V e : (which is done below), we prove the original statement 
as well. 

Let B be a domain such that the surfaces on which Q l _ = are not tangent to the boundary 
B. First, we prove that in the expression 



the integral over {q+} is absolutely convergent (here the integral over {q^} is finite because 
xae > 0, (3 > 0). This becomes obvious (considering conditions (112) and the fact that, in 
type-2 lines, the momentum Q l _ is separated from zero) if the contours of the integration over 
{q 1 !} can be deformed in such a way that the momenta Q\ of the full lines are separated from 
zero by a finite quantity (within the domain V e C\B). In this case, we can repeat the well-known 
Weinberg reasoning [51]. What can prevent deformation is either a "clamping" of the contour 
or the point Q l _ = falling on the integration boundary. 

Let us investigate the first alternative. We divide the domain of integration over q+ into 
sectors such that the momenta of all full lines Q\ have a constant sign within one sector. Let 
us examine one sector. We take a set of full lines whose Q % _ may simultaneously vanish. In the 
vicinity of the point where Q l _ from this set vanish simultaneously, we bend the contours of 
the integration over {q 1 !} such that these contours pass through the points Q l _ = iB % and the 
momenta Q l _ of the type-2 lines do not change. Let B l be such that B l Q l + > for the lines 
from the set (for Q l + from the sector under consideration). It is easy to check that this bending 



2 
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is possible. (Since the contours of integration over q k _ are bent and Q l _ are expressed in terms 
of 5*, one should only check that such b k exist, where the necessary B l are expressed in the 
same way, i.e., that B % obey the conservation laws and flow only along the full lines). With 
this bending, rather small in relation to the deviation and the size of the deviation region, the 
contours do not pass through the poles because, for the denominator of each line from the set 
in question, we have 

(2Q i + Q i _ -M? + 233) -> (2Q\ (Qi + iB*) — Mf + *ae) , Q\B { > 0, (AAA) 

and for the other denominators, the bending takes place in a region separated from the point 
where the corresponding momenta Q l _ are equal to 0. Repeating the reasoning for all sets, we 
can see that there is no contour "clamping". 

The other alternative is excluded by the above condition for B. To make this clear, one 
should introduce such coordinates £ a in the g fc -space that the boundary of the domain B is 
determined by the equation = a = const and then argue as above for the coordinates £ a 
with a > 2. 

After bending the contours, integral (A. 1.3) is absolutely convergent in q+, q k _ if the inte- 
gration in q+ is carried out within the sector under consideration. On pointing out that the 
result, of internal integration in (A. 1.3) does not depend on the bending, we add the integrals 
over all sectors and conclude that (A. 1.3) converges in {q+} absolutely. 

Now let us prove that if B is a quite small, finite vicinity of the point that lies outside 
the sphere B^, then expression (A. 1.3) is equal to zero. We consider the momentum Q l _ of one 
line. Flowing along the diagram, it can ramify or it can merge with other momenta. Clearly, 
two situations are possible: either it flows away completely through external lines, or, probably, 
after long wandering, part of it, Q-, makes a complete loop. The former situation is possible 
only if \Q % _\ < J2 r b-l? where all external momenta leaving the diagram (but not entering it) are 
summed. Obviously, S can be chosen such that for {q 1 !} from B, a line exists whose momentum 
violates this condition. 

The latter situation results in the existence of a loop, where the inequality Q % _ > Q_ holds 
for all momenta of its lines and the positive direction of the momenta is along the loop. Then 
the integral over q\ of the loop in question can be interchanged with the integrals over {q k } 
(because it is absolutely and uniformly convergent for all q^_) and the residue formula can be 
used to perform this integration. Since, for the loop in question, the momenta Q l _ of the lines 
of this loop are separated from zero and are of the same sign, the result is zero. This has a 
simple physical meaning. If we pass to stationary noncovariant perturbation theory, we find 
that only quanta with positive Q- can exist. In this case, external particles with positive p- 
are incoming and those with negative p_ are outgoing. Then, the momentum conservation law 
favors the occurrence of the first situation. 

The entire outside space for B can be composed of the above domains Bl (everything 
converges well at infinity due to the factor exp(— (3J2iQ- ))• Thus, on the left-hand side of 
(A. 1.1), one can substitute the integration domain V e H B L for V e , set the limit in 7 under the 
sign of integration over {q+} because of its absolute convergence, and also set the limit in (3 
under the integration sign because the domain of the integration over {q 1 !} is bounded. Thus, 
we obtain the right-hand side. The statement is proved. 

Statement 2. IfV £ corresponds to the presence of type- 2 lines alone, then, under the same 
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conditions as in Statement 1, the equality 

f(Q\p s ) 



\Q l _ -Mf + iae) 



f(Q l ,P s 



v £ nB L 



? dg -n J (2Q l + e-M2 + ^ae) 



. (A1.5) 



is valid. 

The proof of this statement is analogous to the second part of the proof of Statement 1. 

Appendix 2 

Statement. If conditions (112) are satisfied, the limits inj and (3 in (118) can be interchanged 
(in turn) with the sign of the integral over {a{\ and then with f ■ 
To prove this, we define the vectors 



{q 1 + ,q 1 _,...,q l + ,q l _}^S, . . . ,Q n + ,Ql} = »S + P, 

{yt,yT,---,y^Vn} = Y > 



(A2.1) 



where the vector P is built only from external momenta and \x is an / x n matrix of rank /, 

tAk-l = A*2fe 



• 2l " 2l 1 - 0, nil — /i^ 1 - Next, we introduce the following notation: 



A,: 



7 

— ioti (3 



A = diag{A 1 , . . . , A n }, A — //A/x, 



B = fi'AP - -in*Y, C = -P f AP + iY l P - i ]T a t Mf. 
2 i 

Then it follows from (119) that 



(A2.2) 



Vdet A 



Vi=0 



. (A2.3) 



The function / is a polynomial and we consider each of its terms separately Up to a factor, 
each term has the form -rp— ■ ■ ■ ir—- These derivatives act on C and B. The action on C 
results in the constant factor iN l P, the action on B results in the factor — (l/2)iN t fiA~ 1 B or 
— (l/A)N{nA~ 1 n t N 2 (the latter is the result of the action of two derivatives; N, N 1: and N 2 are 
constant vectors). 

It is necessary to prove the correctness of the following three procedures: (i) setting the 
limit in 7 under the integral sign for fixed (3 > 0; (ii) setting the limit in f3 for 7 = 0; (iii) 
setting the limits in 7 and (3 under the signs of differentiation with respect to Y . In cases (i) 
and (ii), one must obtain the bounds 



\(p(ai,p s ,-f,(3)\ < ip'(a h p s ,/3), 



(A2.4) 
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<p(a i:P s AP)\ <y/(a,,p s ), 



(A.2.5) 



where ip' and ip" are functions integrable (for </?' if f3 > 0) in any finite domain over ctj, with 
ai > 0. Then, for case (i), we have 



i.e., a limit on the integrated function arises, and, thus, the limit in 7 can be put under the 
integral sign. The situation is similar for case (ii). It is evident from (A. 2. 3) that the function 
ip can be singular only if the eigenvalues of matrix A become zero. On finding the lower bound 
of these eigenvalues, one can prove through rather long reasoning that bounds (A. 2. 4), (A.2.5) 
exist if condition (112) is satisfied. 

After the limits in 7 and (3 are put under the integral sign, it is not difficult to interchange 
them with the differentiation with respect to Y. One need do it only for «j > (for each i) 
and, in this case, one can show that the eigenvalues of the matrix A are nonzero and <p is not 
singular. 
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